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Abstract. 

Even if neutrino masses are unknown, we know neutrinos are much lighter than 
the other fermions we know, and we do not have a good explanation for it. In the 
Standard Model of elementary particles neutrinos are exactly massless, although this is 
not insured by any basic principle. Non-zero neutrino masses arise in many extensions 
of the Standard Model. Massive neutrinos and their associated properties, such as the 
Dirac or Majorana character of neutrinos, their mixings, lifetimes and magnetic or 
electric moments, may have very important consequences in astrophysics, cosmology 
and particle physics. Here we explore these consequences and the constraints they 
already impose on neutrino properties, as well as the large body of experimental and 
observational efforts currently devoted to elucidate the mystery of neutrino masses. 
Several hints for non-zero masses in solar and atmospheric neutrinos, that will be 
confirmed or rejected in the near future, make this field of research particularly exciting 
at present. 
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1. Introduction 

Pauli proposed in 1930 the existence of neutrinos as a "desperate way out" to 
explain the continuous spectrum of electrons emitted in /?-decay. This idea was 
considered at the time almost as revolutionary as the alternative explanation, the 
violation of the principle of energy conservation. In this way the neutrino became the 
first particle proposed as the solution to an elementary particle problem, an idea that 
proved to be useful many times later and now commonly used (and even over used). 
Pauli postulated the existence of a new neutral light fermion. Even if we now know 
that neutrinos exists and come in three varieties or "flavours" (as charged leptons do) 
and we have no doubts about their neutrality, we are still debating their masses. 

In fact, fermions come in families or generations, each one a repetition of the 
others, except for their mass. The lighter charged member of each family is heavier 
than the heaviest of the previous one. This is an inter-familial mass hierarchy that has 
no explanation so far, since in the Standard Model (SM) of elementary particles all 
masses are free parameters (see section 3.1). Why the generations are three is another 
mystery. In any case, we know neutrinos come only in the three known varieties, 
unless the additional ones are inert, i.e. do not have weak interactions unlike the 
known neutrinos, or they are heavier than ~ 45 GeV. This results from the 

measurement at the CERN e"'"e~ collider LEP of the width of the weak boson Z° into 
"invisible" particles (neutrinos or other exotic particles that do not interact within 
the LEP detectors), Z° vv. 

At present we only know upper bounds on neutrino masses, obtained in direct 
mass searches (see section 4.1), 

rriy^ < 0(10 eV) m^^ < 160 keV m^^ < 31 MeV . (1.1) 

The bound on my^ is actually uncertain because unknown systematic effects are 
probably responsible for the negative experimentally measured m^^ values. Even if 
there is a combined bound of 5 eV with statistical confidence level of 95%, the above 
mentioned systematic effects bring the bound at about 10 eV. A preliminary result 
lowers the bound on to 29 MeV. These bounds show that neutrinos are much 
lighter than their corresponding charged fermion, me = 0.5 MeV, = 105.6 MeV 
and rrir = 1.77 GeV, and the other members of their respective family. 
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Thus neutrinos introduce an intra-familial hierarchy problem: why neutrinos are 
much hghter than the other members of each family. In the SM this last hierarchy 
is obtained by avoiding producing neutrino masses through the mechanism that gives 
origin to all the other masses, namely the vacuum expectation value of the standard 
scalar Higgs field (see section 3.1). Thus neutrinos are exactly massless in the SM. 
However this is a rather ad-hoc choice. No basic physical principle insures the 
masslessness of neutrinos (as is the case for the masslessness of the photon, insured 
by gauge invariance). 

A large body of experimental and observational efforts are devoted right now to 
elucidate the mystery of neutrino masses. These are direct mass searches, neutrinoless 
double beta decay (PPOv) (or neutrinoless plus a boson (PPOuJ) decay experiments), 
oscillation experiments in reactors and accelerators, solar neutrino observations, 
atmospheric neutrino observations and supernova neutrino observations. 

While direct mass searches do not rely on any other possible neutrino property 
besides its mass, PPOu requires neutrinos to be Majorana particles. Majorana 
proposed in 1937 that neutrinos, contrary to their charged family companions, can be 
their own antiparticle (see section 2). Particle and antiparticle carry opposite charges 
of any conserved lepton number. So if neutrinos are Majorana fermions, lepton number 
is not conserved. This would induce PPOv where lepton number is violated by two 
units. While double beta decay with emission of two neutrinos has been observed, 
PPOu has not. This places an upper bound on an effective Majorana electron neutrino 
mass (see section 4.2) 

(m^J ^ 1 eV . (1.2) 

If neutrinos are massive, there would be mixings between neutrino flavours (see section 
3.1). One of the most striking manifestations of neutrino mixing would be neutrino 
oscillations (see section 6). These appear because of the difference between weak 
interaction neutrino eigenstates and mass eigenstates, if neutrinos are massive. When 
a neutrino is produced it is necessarily in an interaction eigenstate, whose different 
mass eigenstate components propagate differently, according to their mass. Thus, 
after some time the propagating neutrino becomes a linear combination of interaction 
eigenstates, different from the initial one. Actually, because neutrino oscillations are 
an interference effect, they are sensitive to a combination of mixing angles and mass 
square differences. Am?. Oscillations may reveal masses much smaller than those 
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testable through direct searches. Oscillations could explain the so-called solar neutrino 
problem and/or the atmospheric neutrino problem, and are actively being looked for 
at dozens of experiments at present. 

The solar neutrino problem consists of the deficit of observed neutrinos emitted 
from the sun with respect to the theoretically expected amount. It has been with us 
for more than 30 yr, but for a long time it was based on only one experiment, led 
by R. Davis. Several other experiments have now also observed solar neutrinos. One 
of them, Kamiokande, actually became the first neutrino telescope when it succeeded 
in identifying the sun in the sky through its neutrino-image. The results of the four 
experiments that have observed solar neutrinos so far, even if still inconclusive, support 
the existence of a solar neutrino deficit, and suggest that the solution lies in new 
neutrino physics (as opposed to a modification of our standard solar model). This is 
a very exciting field of research and several new experiments that are under way will 
allow us to understand fully the solar neutrino problem within the next few years. 
One of the favourite explanations of this problem involves matter induced oscillations, 
with very small neutrino mass differences, Am^ ~ 10~^ eV^ (see section 7). 

Atmospheric neutrinos produced by cosmic rays hitting the earth atmosphere show 
a deficit of V/j, relative to Ve with respect to the expected ratio. This deficit could be 
explained if part of the oscillated into Vr (oscillations into Ve or a sterile neutrino 
are not allowed, see section 8) within the atmosphere, before reaching the surface of 
the earth. This would require Am^ ~ 10~^ eV^ (see section 8). New long-baseline 
oscillations experiments are being proposed to test experimentally this hypothesis 
(see section 6). The baselines of these experiments are actually amazingly long. For 
example, neutrinos produced at CERN, in Switzerland, could be detected at the Gran 
Sasso Laboratory, in Italy. 

Neutrinos are not only important in particle physics but also in astrophysics 
(section 9) and cosmology (section 5). Neutrino emissions are an important, and 
sometimes dominant, energy loss mechanism in stars. So much so that neutrino 
properties can be constrained on the basis of the unacceptable changes they would 
introduce in the evolution of stars (section 9). Nowhere is the neutrino energy loss 
more striking than in a supernova explosion in which 99% of the energy released goes 
into neutrinos. This is the only case in which the matter densities achieved are so high 
that neutrinos are temporarily trapped. The observation of 19 of the neutrinos emitted 
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by the supernova SN1987A inaugurated neutrino astronomy outside the solar system, 
and brought about an amazing amount of information on neutrinos. For example, a 
bound on the v^, mass was obtained, rriu^ < 23 eV. We would need the observation of 
a supernova in our galaxy (as opposed to SN1987A that happened in one of the Milky 
Way satellite galaxies, the Large Magellanic Cloud) to be able to say something about 
(even measure) the masses of the other neutrinos (if they are larger than about 25 eV 
with the experiments under construction at present). 

Turning to cosmology (section 5), the theory of nucleosynthesis as well as the dark 
matter problem and its related issue of structure formation in the universe, would be 
greatly affected by new neutrino properties. The primordial nucleosynthesis of light 
elements (D + ^He, ^He and ^Li) is one of the fundamental pieces of evidence on which 
the Big Bang model relies. Some neutrino properties, such as a finite mass, lifetime 
or magnetic moment may considerably change the outcome of nucleosynthesis and, in 
fact, this is used at present to constrain neutrino properties. Most of the mass content 
of the universe, 90 to 99%, is in a form of matter that does not emit or absorb light in 
any observable way. This is the dark matter (DM) . Neutrinos of mass in the range of 
a few eV to a few tens of eV could be an important component of the DM. Structure 
formation arguments (namely the formation of galaxies and clusters of galaxies) seem 
to require these light neutrinos not to constitute the bulk of the DM. These neutrinos 
would be hot DM (HDM, i.e. relativistic when galaxies should start forming, namely 
at a temperature of approximately 1 keV) while structure formation prefers the bulk 
of the DM to be cold (CDM). However, CDM does not seem to account by itself for 
all the observations, and an admixture of something else seems to be required. This 
could be some neutrinos as HDM, with mass of a few eV, or heavier unstable neutrinos 
with a tuned combination of masses and lifetimes. 

In this brief introduction we see that many important effects in physics, 
astrophysics, and cosmology depend on neutrino properties, all associated with their 
mass. We have not yet mentioned one of the major motivations that particle physicists 
have to find neutrino masses. This motivation has to do with the triumph and the 
tragedy of the standard model of elementary particles, both joined in the fact that it 
works extremely well. So well that at this point any sign of failure and a consequent 
indication to go beyond this model will be more welcome than a new confirmation. 
This is what any non-zero neutrino mass would be. A potentially rich window towards 
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physics beyond the standard model. Moreover, the expectation of this opening towards 
new physics is intense at present because the solar and atmospheric neutrino problems 
seem to be giving hints of non-zero masses, that will be confirmed or rejected in the 
near future. 

This review is organized in the following manner. We start by explaining the 
different neutrino masses corresponding to different type of neutrinos, Dirac or 
Majorana, in section 2. In section 3 we discuss the main elementary particle models 
for neutrino masses, starting with an explanation of the origin of the masslessness of 
neutrinos in the SM. We present then the main types of proposed extensions of the SM 
that incorporate neutrino masses (plain Dirac masses, left-handed neutrino Majorana 
masses, see-saw models) and their distinctive phenomenological consequences. 

In section 4 we describe the status of direct mass searches and Majorana mass 
searches in neutrinoless double beta decays. In section 5 we go over the many 
cosmological implications of, and constraints on, neutrino properties, mainly masses 
and lifetimes. Sections 6, 7 and 8 review neutrino oscillations, the solar neutrino 
problem and the atmospheric neutrino problem, their implications for neutrino masses 
and the current and future experiments with which these problems will be clarified 
in the near future. We explain oscillations in vacuum in section 6 and oscillations in 
matter in section 7 (as well as other possible solutions to the solar neutrino problem). 
Section 9 describes stars, mainly SN1987A, as a laboratory for neutrino physics and 
summarizes the main bounds they impose on neutrino masses, lifetimes and magnetic 
and electric moments. A few concluding remarks follow. 



2. Types of Possible Neutrino Masses 

Because the electroweak interactions of the Standard Model (SM) violate parity P 
(and charge-conjugation C) maximally, they distinguish fermions of different chirality. 
These are eigenstates of the Dirac matrix 75, with eigenvalue for right-handed 
chirality and —1 for left-handed chirality. Chirality eigenstates are called Weyl spinors. 
Starting from a Dirac spinor ■0, that has four independent complex components, we 
obtain two orthogonal two-component Weyl spinors of definite left and right-handed 
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chirality, ipL and ipn, by means of the projectors Pl and Pr, 



1-75 



1+75 

2 



(2.1) 



2 



The property 7I = 1 insures that Pj, and are, in fact, projectors onto orthogonal 
components, i.e. P| = 1, = 1 and PrPl = PlPr = 0. The decomposition. 



is invariant under the homogeneous Lorentz group (i.e. the continuous transformations 
of the Lorentz group, rotations and boosts), because "(pL and transform as separate 
irreducible representations. However, a discrete parity (P) transformation transforms 
left and right-handed Weyl spinors into each other. This is why a P-invariant theory 
must treat both chiralities identically (equal couplings) and conversely, when only one 
chirality is present, the violation of parity is maximal. The experimentally determined 
V — A (vector minus axial vector) coupling in charged current weak interactions means 
that only left-handed fermions interact with the gauge bosons of the group SUl{2) 
and P-violation is maximal. For example, the leptonic charged current term in the 
Lagrangian of the SM is 



Here g is the weak coupling constant, £ = e, are the charged leptons of the known 
generations, ui are the associated neutrinos and W~ is the charged heavy gauge boson. 
Notice that 



because 75 anticommutes with all the other Dirac matrices 7^ = 70)7i)72)73) i-e. 
757m = -7/u75, and, thus ipR = ipPL, ipL = ipPn- 

A Weyl spinor can describe only free massless particles, for which chirality coincides 
with helicity (the spin projection in the direction of motion), that is a good quantum 
number. This is so because mass terms mix chiralities. Thus, if a fermion is massive, 
the helicity eigenstates, the eigenstates of a freely propagating fermion, do not coincide 
with the chirality eigenstates, the eigenstates of weak interactions. A massive fermion 
can be either a Dirac or a Majorana particle, that are described respectively by a 
Dirac and a Majorana spinor. 



(2.2) 




(2.3) 




(2.4) 
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A Dirac spinor, such as the one describing an electron, has four independent 
complex components, corresponding to particle and antiparticle, both of left- and 
right-handed helicity. It was actually Dirac's equation for the electron that led to 
the concept of particles and antiparticles and to the definition of charge-conjugation. 
Antiparticle spinors are obtained from particle spinors through charge-conjugation, 
under which il) — > rycV''^; with r]c a phase and 

r = C^^ = C^or- (2.5) 

Here C is a unitary matrix that is defined by the property C~-^7^C = — 7^- If we 
now apply Pl or to the Dirac spinor and its charge-conjugate spinor, we see that 
under the physical charge conjugation operation, that should preserve chirality* (and 
ignoring the phase r]c hereafter), t/j^ — ^ ('*/''^)l and ifjfi {iI^'^)r (see, e.g., Langacker 
1981). Hence, under this operation, a Dirac particle in a given state of momentum and 
helicity is changed into its antiparticle in the same state of momentum and helicity. 

Notice that one does not obtain the charge-conjugate of the Weyl spinor 
components of a Dirac spinor by applying the conjugation operation (2.5) directly 

to them. In fact, the operations of conjugation (2.5) and projection onto chirality 

J' 

components do not commute, (V'l)'^ = CipL is a i?-handed spinor, it is the right 
component of ip'' in (2.5), (iPl)^ = {'4'^)r- Hence, in this operation the chirality is 
reversed, while, as discussed above, under charge conjugation i/jl = (V'i?)'^ 

(thus, an interaction containing only L-handed fields violates charge-conjugation 
maximally). Many authors then understandably avoid calling {i^lY ^^e charge- 
conjugate of ipL. They call it simply "conjugate" (see e.g. Mohapatra and Pal 
1991), or even loosely "CP-partner" (see e.g. Langacker 1981 and 1992), because 
of the additional change of chirality, even if it is not the CP-conjugate**. We will 
then call c-conjugation the operation defined in (2.5) and we call '4>'^ the c-conjugate 
of ij) for any kind of field ip, Dirac or Weyl (or later Majorana). The operation c- 
conjugation coincides, up to a phase, with charge-conjugation C for a Dirac field (and 
for a Majorana field) but not for a Weyl spinor component of a Dirac field. It is 

* It should preserve helicity, that is the actual physical property of a particle, but thinking of the 

limit in which chirality and helicity coincide, it is obvious that chirality must also be preserved. 

** Actually CP conjugation is a different operation, involving an additional multiplication by 70 and 

a change of x into —x, ■4>^^{t,x) = »7cp7oV'"^(*j where t]cp is a phase. 
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also trivial to check that using (V'l)'^ = (V''^)fl {'4'rT — ot using the C- 

transformation ipL ^ i'^^)L and ipn — > {iP'^)r yields the same result when applied to 
V'L + ipR (we only exchange which component of if^'^ we call the conjugate of which 
component of 'i/')- 

In the language of relativistic quantum field theory, a Weyl field, say V'l, can 

annihilate a left-handed (L) particle or create a right-handed (R) antiparticle, while 
T 

= Ci/^L -I can annihilate a R-antiparticle or create a L-particle. For the roles 
of right and left are exchanged. Notice that V'l and {'4>'^)r are not independent fields, 
they represent the same two degrees of freedom. The same is valid for ij^R and {iP^)l- 
For Dirac fermions ipL and ipj^ represent different degrees of freedom but for Majorana 
fermions they coincide. Since the fields ipL and {iP^)l must have opposite values of all 
additive quantum numbers, Majorana fermions cannot carry any conserved charge. 

Electrons e and positrons e'^ are clearly different particles, since they have opposite 
electric charges. However, wc have examples of neutral particles, such as the iv^, that 
coincide with their antiparticlcs. Majorana (Majorana 1937) first proposed that a 
neutral fermion could have this property. While a Dirac fermion is always diff'crent 
from its antiparticle, a Majorana fermion is such that (ip^)^ = ip^ , or actually, 

since a phase, e*® (called sometimes Majorana creation phase factor, Kayser 1984), 
can always be incorporated in the definition of the fermion i/j^ . Instead of 
four independent components, a Majorana spinor has only two, since particle and 
antiparticle coincide. 

Is u'' different or identical to u? 

Only neutrinos of left-handed chirality and antineutrinos of right -handed chirality 
(present in the hermitian conjugated terms of (2.3)), and {i''^)r, interact in the 
SM. If they exist, the components of opposite chirality, ur and (z^°)l, are "inert" 
since they do not participate in weak interactions. If neutrinos are massless they are 
described by Weyl spinors, chirality and helicity coincide and we can never produce 
neutrinos and antineutrinos of the same helicity in weak processes. "Neutrino" and 
"antineutrino" are just names at this point, for states that interact differently due to 
their different helicity. For massive neutrinos we can change the helicity by a boost 
that inverts the momentum. We can then produce neutrinos and antineutrinos of the 
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same helicity and in principle compare them. If they interact differently, neutrinos 
are Dirac particles, if they are identical, neutrinos are Majorana particles (for a more 
detailed explanation, see e.g. Kayser 1985). Although Majorana and Dirac neutrinos 
have different properties, the differences vanish with the neutrino mass. The only 
feasible experiment that may actually prove the Majorana nature of neutrinos is the 
neutrinoless double beta decay (see section 4.2). 
A Dirac mass term, of the form 

-£^ass = m^^^ = m^{ii;L + ipR){ipL + ^r) = m^i^L^R + ^r^^l) (2.7) 

mixes two different Weyl spinors, of opposite chirality. This is the type of mass 
generated in the SM. Thus, the exclusion in the SM of the vr (and (j^'^)l) components 
insures ad-hoc that neutrinos are massless. A Majorana mass term can be written 
with only one Weyl spinor, say ipL-, and its c-conjugate {i^lY = C{'ipL) ^ = {'4'^)r, 

= \m^W^¥Lr){i^L + ii^LY) . (2.8) 

The mass eigenstate field ip'^ , 

= + (V'l)' = V'l + mR ■ (2.9) 

has the canonical mass term ^m^ip^ip^ and is a Majorana field with creation phase 
e^^ = +1, since (ip^Y = (notice that ((iPlYY = V'l)- The field 

^^ = ^L- {^lY (2.10) 

is also a Majorana field, as defined in (2.6), with the phase e*^ = —1, i.e. {ip^Y = 
—ip^. Clearly ip^ and ip^ describe the same degrees of freedom. It is straightforward 
to check that if we were to use tp^ instead of tp^ in (2.8), the sign of the mass term 
(2.8) would be reversed. 

Notice that a Majorana mass term carries two units of fermion number since 
{iPlY same fermion number as ipL. A Dirac mass term has instead zero 

fermion number, since ipj^ and ipL have opposite fermion numbers. Thus, Majorana 
mass terms are forbidden if the fermion number, the lepton number in the case of 
neutrinos, is conserved. 
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mass — "T?^ 



M 



When both Majorana and Dirac masses are present, in general the massive 
fermions are Majorana particles, and fermion number is not conserved. The most 
general mass matrix using the Weyl spinors vl and is: 

->Cmass =mD (WJ^R + J^I'l) + 



This expression is also valid for an arbitrary number of flavours, just taking uj^ = 
{^eL, ^fiL, ■ ■ •) and similarly for vji. Due to the anticommutation properties of the 
fermion fields, the mass matrix in (2.12) turns out to be symmetric (see e.g. Bilenky 
and Petcov 1987). We recall that a symmetric matrix M can generally be brought to 
a diagonal matrix with positive entries, by means of a unitary transformation U, 
i.e. U'^MU = Ma. 

With one neutrino species, , and are numbers (thus, m^=m£)). Let 
us consider this case first. When = and = 0, the mass eigenstate is a Dirac 
fermion (vl + J^r) with mass . 

In this case the diagonalization of the matrix of (2.12) yields two chiral eigenstates 
ipL± = (i^L± (i^rY) / with equal and opposite mass eigenvalues, mo and —ijid- In 
terms of the associated Majorana fields ip± = ipL± + {''Pl±Y, that are self-conjugate by 
construction, the mass term reads —2Cmass = rnoi^-^-i^'l. ~ ''^ni^-i^-- The negative 
sign in the mass can however be absorbed into a redefinition of the corresponding 
Majorana eigenstates, taking them to be t/ji = tj^j^ and '02 = V'l- ~ [ipL-Y, and this 
last will consequently transform under conjugation as the field i/j^ discussed above. 
Notice that the original sign of the mass, call it A^, becomes now the Majorana 
creation phase (the phase in (2.6)) of the mass eigenstate field with positive mass. We 
also note that Majorana particles, unlike Dirac ones, have a well defined intrinsic CP 
parity, fjcp*- This is so because CP transforms a particle field into the antiparticle 
field and while for Dirac neutrinos there is the freedom to redefine the phase of the 

* VCP — e^^VCP — i^) where rjcp is the phase of the CP-transformation (see footnote in previous 
page) and e*^ is the Majorana creation phase. 



(2.11) 



In matricial form it becomes (using vlvh = {I'RYii'LY) 




(2.12) 



13 



antiparticle field to absorb the CP phase, this is not possible in the Majorana case. 
It is always possible to chose the neutrino fields, as we did in our example, so that 
fjcp = Xmi (see e.g. Kayser 1988 and 1984 and Bilenky and Petcov 1987). 

The linear combination of degenerate eigenstates v = {ipi + 1^2) /V2 is a Dirac 
neutrino, u — {ul + i^r), and u'^ = (-01 — 1^2) /V^)- Thus, a (four-component) Dirac 
field consists of two (two-component) degenerate Majorana fields with opposite values 
of (or, in terms of particle states, a Dirac state can be viewed as the sum of two 
degenerate Majorana states with opposite CP-parities). The phase is irrelevant 
for freely propagating neutrinos, but it will appear in the charged currents involving 
the {i'lY field, affecting for instance the neutrinoless double beta decay (see section 
4.2). 

These phases appear generically when dealing with Majorana neutrinos. The 
reason is the following. When CP is conserved, arbitrary phases in the original fields 
can be chosen so that the mass matrix is real. Because it is also symmetric, and 
hence hermitian, one is tempted to diagonalize it by means of an orthogonal (not 
just unitary) transformation O. The price to pay is, then, the possible appearance of 
negative eigenvalues for self-conjugate fields or antiself-conjugate fields, if one wants 
positive masses (as in the simple example of two degenerate eigenstates presented 
above). An alternative procedure would be to insist in having self-conjugate fields 
with positive masses, at the expense of having a non-real transformation matrix. This 
is obtained by just replacing the if^nL = J^^ni^iL mass eigenstates with negative 
masses (here n labels mass eigenstates and / interaction eigenstates) by the states 
i^nL = ^i^nL = Yl,{iOni)nL- The Majorana states ^i'^ = •i/'nL + {'^'ulY ^iU be self- 
conjugate and have positive masses, but their mixing matrix (and their coupling to 
the charged gauge bosons) will involve the non-real matrix Uni = iOni- These different 
approaches are of course generalizable to the CP violating case. 

Going back to (2.12), when either or or both are non zero, the mass 
eigenvectors are two Majorana spinors with different mass. When the Majorana 
masses are small compared to the Dirac masses, the resulting eigenstate is a "pseudo-" 
Dirac neutrino (Wolfenstein 1981), in which the two Majorana fields are not exactly 
degenerate but almost so. 

In order to incorporate the three known generations, we take each of the spinors 
to be a vector in generation space, = (j^blj ^ij-li ^tl)j ^^c. and , and to 
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be matrices. If there are as many right-handed neutrino species as left-handed ones, 
the sub- matrices in (2.12) are 3x3. However, it may be that uji are entirely absent, 
so only exists, or that the number of vr is different from three. In the case of 
several generations there is a special type of Dirac neutrino that may appear as a mass 
eigenstate, the ZKM type (Zeldovich 1952, Konopinski and Mahmoud 1953). Usually 
a Dirac neutrino left-handed component is "active" (it interacts weakly) while its 
right-handed component is "sterile" (or inert). Both components of a ZKM neutrino 
are active. The right handed component is the antiparticle of the left-handed neutrino 
of a different generation, e.g. u = v^l + {i^ixlY- The ZKM neutrino appears when 
a linear combination of flavour lepton numbers different that the usual total lepton 
number {L = -\- -\- Lj. leading to a usual Dirac neutrino) is conserved, in the 
example it is — L^. This type of neutrino requires a complicated Higgs structure 
not necessary for usual Dirac neutrinos. 

A Dirac neutrino has, in general, a magnetic dipole moment, and may also have 
an electric dipole moment (that violates CP). So, in principle, a torque exerted 
by an external magnetic or electric field, B or E, can change its chirality (since 
dipole moment interactions couple Weyl spinors of opposite chiralities) . Majorana 
neutrinos can only have "transition" dipole moments, (Schechter and Valle 1981 and 
1982, Kayser 1982, Nieves 1982), which in a 5 or ^ field cause the neutrino not only to 
change chirality, but also to change flavour (i.e. a neutrino of one generation is changed 
into another different neutrino of a different generation), for example z/^l — > (^'e)ji- 
Transition moments can exist for Dirac neutrinos too, and they imply also the radiative 
decay of the heavier neutrino into the lighter one and a photon. 



3. Main Elementary Particle Models for Neutrino masses. 
3.1. The Standard Model (SM) 

In the SM all charged fermions, quarks and leptons, get Dirac masses through the 
Higgs mechanism that breaks spontaneously the SUl{2) x Uy{1) electroweak gauge 
symmetry into ?7e.m.(l) of electromagnetism, giving mass to three of the four gauge 

15 



vector bosons that mediate electroweak interaction, the W'^,W~, and while only 
the photon stays massless. The absence of the uji prevents the appearance of a Dirac 
mass for the neutrinos. The accidental global symmetry Ub-l{^) corresponding to the 
conservation of B — L, Baryon number minus Lepton number, prevents the appearance 
of Majorana masses for neutrinos. Thus neutrinos are massless in the SM. 

All masses in the SM, are proportional to the vacuum expectation value (VEV) of 
a complex SUl{2) doublet scalar field (f) and its conjugate (f) = ia'^cf)* (where cr^ is the 
second Pauli matrix), that is also a doublet (because the representations of SU{2) are 
real) , 

A potential energy for is introduced, that has its minimum at 

Since the component 0° that has non-zero VEV is neutral, ?7e.m.(l) is preserved as a 
good gauge symmetry. Given the masses of the W and Z particles, one obtains 

V ~ 250 GeV . (3.3) 

The matter fields are 15 Weyl spinors, repeated twice at larger masses yielding three 
generations or families (see the table 1). The quarks q and leptons £ of the different 
generations are distinguished by their "fiavour" . There are three flavours of charged 
leptons, the electron e, the muon and the r, and their corresponding neutrinos 
i/eji^Hji^T- Quarks carry also "colours", the quantum numbers of strong interactions. 
There are three colours for each quark flavour, since quarks belong to the fundamental 
representation of SUc{^)- In the table 1 the parentheses show the representations of 
SUc{^) and SUl{2) (both indicated by their dimensionality), and the value of the 
weak hypercharge F, for each set of flelds. 

The left handed fermions /l {ql and £l) have weak isospin Tl = 1/2 (T3 is +1/2 
and —1/2 for the upper and lower components). The right handed fermions fa, have 
Tl = 0, (thus TsL = 0). The charge Q of each field is given by Q = T3L + Y. The 
charges of the particles in each generation sum to zero, S/Q = 0. This insures that 
anomalies cancel within each family and, thus, the SM symmetry is preserved at the 
quantum level. 
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Table 1. Matter fields in the Standard Model. Here, i = 1, 2, 3 is the colour index. 



Quarks 

Ql: (3,2,1/6) 

ur: (3,1,2/3) 

URi , CRi , tRi 



e 



Leptons 

(1,2,-1/2) 



T 



dR-. (3,1,-1/3) 

d-R^ , Sfli , hRi 



en - (1,1,-1) 

■, [i-R Tr 



The SUl{2) x ^/^^(l) assignments of and Jr make it impossible to have gauge 
invariant (i.e. gauge singlet) mass terms, i.e. 



mfhfR , 



(3.4) 



in the Lagrangian. Notice that these Dirac mass terms have a weak isospin = 1/2, 
thus they need to be coupled to a weak doublet to form a singlet. In fact, the fermion 
masses come from gauge invariant Yukawa type couplings with the doublet scalar 
0= (1,2,-1/2), 

—-^Yukawa = ^ ] ^{^u 

a, 13 

(3.5) 

where a, [3 = 1,2,3 are generation indices and {Xf)a/3 are coupling constants. 
Replacing by its VEV (3.2), one obtains the fermion mass matrices 



(rn/)a/3 = (A/)a/3;^ • 



(3.6) 



In the SM, the couplings {Xf)ai3 are arbitrary complex numbers, thus all fermion 
masses are free parameters. As we mentioned before, it is only the absence of the 
right-handed neutrinos V(,Ri i'hr-, i'tR-, that prevents the appearance of Dirac neutrino 
masses through Yukawa couplings with the Higgs field. 
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Therefore, in the SM neutrinos are massless because there are no right-handed 
neutrinos and there are no other Higgs bosons besides the doublet (that could 
generate Majorana masses). We could entertain the idea that even in this case 
neutrinos could get masses, Majorana masses, due to higher order corrections. This is 
not possible in the SM due to the accidental conservation of B — L. Actually B and 
the separate flavour lepton numbers Le, and Lj. are accidental symmetries, but 
they are all anomalous symmetries. Only the combination B — L is anomaly-free and 
hence it is a good global symmetry of the SM. This symmetry prevents the appearance 
of Majorana mass terms also beyond the tree level, because they would violate L and, 
therefore, B — L, by two units. 

We can see now that models of non-zero neutrino masses necessarily add to the 
SM either fermions, typically ur, or bosons, or both, and when the masses are of the 
Majorana type, they introduce a violation of B — L, either explicit or spontaneous. We 
can also see here that if neutrinos are massive, the leptonic charged currents in terms 
of mass eigenstates include a mixing matrix Ki equivalent to the Cabibbo-Kobayashi- 
Maskawa (Cabibbo 1963, Kobayashi and Maskawa 1973) mixing matrix Kq present in 
the weak charged currents of quarks. The mass matrices (3.6) must be diagonalized, 
thus quarks and charged leptons are expressed in the basis of mass eigenstates 

ul U'^ul , dL W^dL , bl U^L^L ■ (3.7) 

Because neutrinos are massless in the SM, C/| is absorbed in the redefinition of the 
neutrinos so that one obtains the coupling in (2.3) and (3.8) where Kg = 1, while 
the matrix Kg = Uf^U'^ remains in the charged current of quarks (we ignore colour 
indices) , 

(d\ /e\ 



= {uct)L l^Kq 



Neutrino mass terms in the Lagrangian are always written in the basis of "current 
eigenstates" (or interaction or flavour eigenstates), that are those that diagonalize the 
lepton weak charged current with the charged lepton mass eigenstates (as in (2.3) and 
(3.8)), namely VeT^HT^r- Moreover charged lepton current eigenstates are chosen to 
coincide with the charged mass eigenstates so that C/| = 1. Then, when neutrinos are 
massive, current neutrino eigenstates differ from mass eigenstates by a mixing matrix 
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(3.8) 



Ki = Ul. Thus in (3.8) 

(^current)T ^ (^^ ^ (l^Tn^ K*, = {v^ ...)r^K*, . (3.9) 



where the dots indicate that there could be more than three neutrino mass eigenstates 
(either heavier than ~ 45 GeV or mainly consisting of inert neutrinos). 

3.2. Plain Dirac Masses 

This possibility corresponds to having only mo non-zero in the general mass matrix 
of (2.12). Right-handed neutrinos can be easily added to the SM, as singlets of the 
gauge group SUc{S) x SU{2) x t/y (1), ur : (1, 1, 0). Since they are inert, they do not 
contribute to anomalies and the renormalizability of the model is not affected. We 
could then add neutrino Yukawa couplings to (3.5) (similar to the up-quark terms), 



that would yield Dirac neutrino masses as (3.6), {m,y)af3 = {^v)ai3 v/y/2. This is the 
simplest way to get non-zero neutrino masses. However, this mechanism provides no 
explanation for the smallness of neutrino masses with respect to the other fermions of 
the same generation. All we could say is that {X^)ai3-i being arbitrary as all the other 
Yukawa couplings, happen to be much smaller than the others. Other mechanisms 
provide some insight into this question. 

3.3. Models Without ur 

Counting only with left-handed Weyl spinors, ul, we can only have the Majorana mass 
terms of mf^ in (2.12). These terms have weak isospin Tl = 1. Thus, renormalizable 
interactions at tree level give origin to these masses (through a Higgs Mechanism) only 
if a triplet Higgs field A = (A+''", A+, A°) is added to the SM, with Yukawa couplings 



Here a = (cri, cr2j C3) are the Pauli matrices, £aL = {j^aL, ^ah) are lepton doublets, 
thus {EollY = ii^aLy, -{i^aLY), sincc (^ql)" IS also the SU(2) conjugate of iaL (as 




(3.10) 



(3.11) 
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4> is the conjugate of in (3.1)), a,/? are generation indices, A = (Ai,A2,A3) and 
A° = (Ai +zA2)/V2, A+ = A3 and A++ = (Ai -zA2)/v/2. These Yukawa couphngs 
contain terms 

9c.(3 V2 {p^lYpul^^ + • • • (3.12) 
that yield Majorana masses proportional to the VEV (A°) = vt/^/^ , 

{jn^)al3 = 9apVT ■ (3.13) 

Here the smallness of the neutrino mass depends on choosing ad- hoc a small vt, 
vt -C v. 

When Majorana masses are present in a model there is a choice between explicit 
or spontaneous L violation. Moreover L (actually we should speak of S — L, because 
L is anomalous) may be a gauge or a global symmetry. If L is a global symmetry, 
its spontaneous breaking would generate a Nambu-Goldstone boson, that is called 
Majoron. 

The lepton number of A is La = —2, as defined by (3.11). Therefore, if L is 
conserved in the Lagrangian, (A'^) 7^ violates L-conservation spontaneously and a 
Goldstone boson results. This is the basic mechanism of the triplet Majoron model 
(Gelmini and Roncadelli 1981), now experimentally rejected (see below). 

Given the matter content of the SM (table 1), with no ur, there are only two other 
scalars, besides the standard doublet and the triplet A, that can have /S't/L(2) xt/y (1) 
invariant Yukawa couplings. These are, a singly charged singlet scalar field h'^, first 
proposed by Zee (Zee 1980), that could couple to lepton doublets 

faf3{iaLY£/3Lh+ = /a^(fc}^i^/3L - (^^e^^)/i+ , (3.14) 

and a doubly charged singlet scalar , that could couple to lepton singlets 

Y {ei3R)h++ . (3.15) 

Here, the couplings foi/s and /^^ must be antisymmetric in a and /3 due to Fermi 
statistics {a, (3 denote generations). Thus /i"*" and h'^'^ couple only to two leptons 
of different families. Neither h'^ nor h~^~^ can have a non-zero VEV (otherwise 
electromagnetism would be a spontaneously broken symmetry). Thus L violation 
has to be introduced in the Higgs potential, by adding at least a second doublet Higgs 
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field, and Majorana masses for neutrinos are obtained through radiative corrections 
(Zee 1980). While the Zee model, based on does not include the h'^'^ field, models 
with do require /i+ (Babu 1988). In both models neutrino masses are small 
because they are generated through radiative corrections, and a distinctive feature is 
the presence of a neutrino much lighter than the other two. 

Some L-violating physics beyond the SM, at a large mass scale A, may induce 
fL Majorana masses through effective (non-renormalizable) couplings. The leading 
non-renormalizable L-violating operator that can be written with the fields of the SM 
is of dimension four (Barbieri et al 1980) 

>^eff = Cap 

When the doublet gets a VEV (3.2), (3.16) gives Majorana masses 

{mf!U - ^ . (3.17) 

For A ^ V one naturally obtains small neutrino masses. As we will see below, in 
see-saw models A is the vr Majorana mass scale, A ~ rn^ . 

We may not need new physics to have terms as (3.16). We know that gravity 
is not incorporated into the SM and it may not respect any global symmetry (for 
example, global charges disappear when falling into a black hole). Quantum gravity 
effects may violate L. In this case A ~ Mpianck — 10"*^^ GeV and the resulting masses 
are mf ~ ClO"^ eV. 



(3.16) 



3.4. See-Saw Models 

The "see-saw" mechanism (Yanagida 1979, Gell-Mann et al 1979) consists of making 
one particle light at the expense of making another heavy. It assumes a hierarchy in 
the values of the different elements of the mass matrix in (2.12), namely 

= M :^ mo = jji , (3.18) 

with |U either zero or negligible. 

Consider first the case of just one generation. Then, the mass matrix has one heavy 
eigenvector, mainly consisting of the inert zv^, N ~ [i'R+{i'RY] + {mD/M)[vL + {i'LY], 
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and one light eigenvector, mainly consisting of the active Weyl spinor v ~ 
{^L - {^lY\ + ("7,D/M)[zy_R - {vYi)% with masses 

O 

niN ^M, ~ < mo . (3.19) 

Thus in the see-saw mechanism the larger M the lighter is the light neutrino, explaining 
the intra- familial hierarchy m,^ <^ tud, while the Dirac neutrino masses are naturally 
similar to the Dirac masses of the other fermions of the same generation. 
If is not negligible, the light mass eigenvalue becomes* 



M 



(3.20) 



Thus, unless /j, < mj^/M one loses the natural explanation of the smallness of rriu. 

In the case of three generations, each fermion in (2.12) becomes a vector and 
M, niD and n become matrices. Heavy (N) and light {u) eigenvectors are found 
diagonalizing the matrices 

rriN = M, my — fji — mDM~^rri^ . (3-21) 

Here, for the heavy masses we kept only the leading term. In the usual see-saw 
models there are as many right- as left-handed neutrinos, M, mo and ^ are 3x3 
matrices (with fi negligible or zero (3.8)) and the symmetric 6x6 mass matrix has 
six eigenvectors, (usually) three light and three heavy, with masses 



"^■Dl ^D2 

Ml ' Mo ' M^ 



, Mi,M2,M3 , (3.22) 



where Mi, i=l,2,3, are the eigenvalues of the matrix M. 

There are many different versions of this mechanism. For example, the Dirac 
masses moi in the eigenvalues (3.22) could be similar to the up-quark masses (usual 
in Grand Unified models), or to the charged lepton masses of the same generation (or 
even something diff'erent). The neutrino mass hierarchy depends on this choice since 
mt/mc ^ 100 while mr/m^ ~ 17. In a "quadratic see-saw", the three heavy masses 

* When this eigenvalue is negative before taking the modulus, the minus sign is incorporated in the 
definition of the light eigenstate, that becomes anti-selfconjugate (see section2), as is the case of v 
above. 
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Mi are similar, i.e. Mj ~ M for i = 1,2,3. Consequently, the hierarchy of light neutrino 
masses is that of m|,^, i.e. m^j^ : rrii,^ : rui,^ ~ : : m| (or : rn^ : m^). In 
a "linear see-saw", the hierarchy of the heavy masses Mj coincides with that oi niDi, 
Ml : M2 : M3 ~ mu rric : mt (for example, Mi generated through loops in a SO(IO) 
model, Witten 1990). As a consequence, also the ratio of light neutrino masses is 
linear in m/j. Vastly different scales are possible for M, depending on the model, 
from 1 TeV or lower (for example, in left-right symmetric models), to 10^^ GeV or 
higher (in GUT's). Only in the simplest GUT models the leptonic mixings coincide 
with those of quarks, K,^ = Kq (up to renormalization effects) in (3.8) and (3.9), but 
one usually also obtains incorrect relations of mass ratios (namely me/m^ = rrid/ms). 
Otherwise the mixing angles are entirely model dependent. 

There are many other less usual possibilities for see-saw mechanisms. For example, 
in the "incomplete see-saw" (Johnson et ai 1986, 1987, Glashow 1991), a 3x3 matrix M 
but of rank 2, i.e. with one zero eigenvalue, yields a Dirac neutrino mass eigenstate of 
mass ~ m^), besides two heavy and two light neutrinos, of masses ~ M and ~ m\,/M 
respectively, as before. 

3.5. Majoron Models 

When a global continuous symmetry is spontaneously broken, a zero mass boson, 
called Nambu-Goldstone boson, appears for every broken generator. One may wonder 
if such massless bosons would lead to new gravitational or Coulomb- like potentials, 
against which there are very stringent upper limits. This is not the case. Goldstone 
bosons generate only spin-dependent non-relativistic long-range potentials for 
which bounds are much less restrictive (Gelmini et al 1983). 

Majorons are the Goldstone bosons associated with the spontaneous breaking of a 
leptonic global symmetry, at a scale V, that usually also generates Majorana neutrino 
masses (thus the name Majoron), rrii, ~ g^ffV. The effective coupling g^ff is also 
roughly the coupling of neutrinos with the Majoron, J. These couplings provide 
the most important phenomenological consequences characteristic of these models. 
In general, given a certain range of masses mj^, for small values of V (typically 
V < V 100 GeV), the couplings g^ff are large, neutrinos may decay {vh — > viJ) 
and/or annihilate {vv — > J J) or interconvert {vv — > v'v') very fast (so neutrinos in the 
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early Universe could not remain to be the Dark Matter). Besides, the Majoron could 
be emitted in neutrinoless double beta decays at a level that could be observable soon 
and Majorons may play a role in energy loss mechanisms in stars or in the collapse 
of a supernova. If instead, given a certain range of masses mj^, V is large (typically 
V > V 100 GeV), then g^ff is relatively small and the only important consequence of 
Majorons is, usually, that neutrinos may decay with much longer lifetimes than before, 
shorter or longer than the lifetime of the universe (so that, for example, neutrinos could 
be a relevant part of the Dark Matter) . 

Majoron models are almost unique in allowing neutrino masses in the range 
between ~ 100 eV and a few GeV, otherwise forbidden due to cosmological arguments, 
through neutrino decays or annihilations into invisible modes in the early universe. 
Another important consequence of Majoron models is that the standard Higgs particle 
may be "invisible" itself through decaying dominantly into the invisible channel 
H J J (Schrock and Suzuki 1892, see for example Lopez-Fernandez et al 1993 
or Brahmachari ct al 1993 and references therein). 

The main properties of these models are determined by the weak isospin of the 
Majoron. The singlet Majoron model (Chikashige, Mohapatra and Peccei 1980, 1981), 
adds to the SM right-handed neutrinos and a singlet Higgs field a coupled to them, 

hal3{l^aRYt^f3RCr + h.C. . (3.23) 

The VEV (cr) ^ 0, generates Majorana masses for the vr, M^p ^ ha/3{cr). Light 
neutrino eigenstates Va result from the see-saw mechanism ((2.12) with = 0), 
if (cr) is large enough, (cr) ^ v. Because L is chosen to be an exact symmetry of 
the Lagrangian and (3.23) defines a non-zero lepton number for cr, Lcr = —2, L is 
spontaneously broken. The Majoron in this model is J = \/2Im{a), thus 

a^^[V + p + iJ] . (3.24) 

Thus, the Majoron couple to the light neutrinos only through the small admixture of 
ur, of order {uid/M), that they contain. Consequently, the decay rate of a heavier 
light neutrino i^h into a lighter one vi, is 

r(i/, - J) ^ (^)^m., . (3.25) 

Actually, in the simplest form of the singlet model presented here, the leading 
(m£)/M)^ terms in the amplitude can be rotated away (Schechter et al 1982), so 
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that r ~ {mD/M)^m,^f^. Thus v hfetimes are large with respect to the hfetime of the 
universe. This is not true in more comphcated versions of the model (Gelmini and 
Valle 1984, Jungman and Luty 1991, Babu 1991, Glashow 1991). 

Besides its effects on neutrinos, a singlet Majoron is practically invisible. It has 
very small couplings to charged fermions, only through neutrino loops, and, because 
the field cr is a gauge singlet, J is not coupled to the boson, thus it does not 
contribute to its invisible decay modes. 

Non-singlet Majorons have been rejected by the LEP bound on the effective 
number of neutrinos, N„ < 2.983 ± 0.025 (Review of Particle Properties 1994). The 
reason is that in these models there are additional contributions to the invisible width 
of the Z°-boson that would count as extra light neutrinos. The new decay mode is 
Z pj, where p is a light boson associated to the Majoron (usually p is the real 
part of the same combination of fields of which J is the imaginary part, as in (3.24)). 
Because non-singlet Majoron models are less "invisible" than singlet Majorons, the 
scale of L- violation V is phenomenologically required to be ^ v ~ 100 GeV, and 
in this case rup <S V, and p can be emitted in a Z^ decay. 

The models rejected are the triplet Majoron (Gelmini and Roncadelli 1981, Georgi 
et al 1981), the doublet Majoron (Bertolini and Santamaria 1988) and supersymmetric 
Majoron models where it is the left-handed scalar neutrino ul (the supersymmetric 
partner of u^) VEV that violates L spontaneously, V = (ul) (Aulakh and Mohapatra 
1983, Ross and Valle 1985). In these three models the additional Z^ width, r{Z^ 
pJ), equals 2, 0.5 and 0.5 respectively of its partial width into a light neutrino 
species, r(Z'° ^a^a)j while only a few percent is still allowed by the LEP results. 
Therefore viable Majorons are singlet (coupled or not to uji) or mostly singlet (even if 
mixed with non-singlets). Examples are the original CMP singlet model (Chikashige, 
Mohapatra and Peccei 1980, 1981), the singlet-triplet or "invisible triplet" (Choi et al 
1989, Choi and Santamaria 1991, D'Ambrosio and Gelmini 1987) and supersymmetric 
models where the right-handed sneutrino (an electroweak singlet as its supersymmetric 
partner lyji) VEV breaks L spontaneously, V = {vr) (Masiero and Valle 1990). 

The "invisible" triplet model preserves many of the characteristics of the original 
triplet Majoron. This is a variation of the triplet model presented in section 3.2, 
which yields a Majoron if L is only spontaneously violated by (A*^) ~ 7^ 0. The 
triplet Majoron is mainly J ~ /mA°, thus it couples at tree level to SM neutrinos. In 
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the "invisible" triplet a singlet a (not coupled to fermions) is added, whose VEV 
(cr) ~ 7^ also breaks L. If vg > vt the Majoron is now mostly singlet, 
J ~ Im{a + {vt/vs)A^)-i and p is the real part of the same combination. Thus 
now V{Z — > pj) ~ (vt/vs)'^ is reduced to an acceptable level. 

Singlet models tend to have large V (and thus small Qes)- However, several singlet 
and mixed models have been proposed (Berezhiani et al 1992, Burgess and Cline 1993, 
1994, Carone 1993) in which ^feff is large (and V small), mainly motivated by the 
possibility of finding Majoron emission in neutrinoless double beta decay, i.e. having 
{gef{)i'eVeJ ^ot much smaller than its experimental upper bound of 0.7 x 10~^ (Moe 
1994, see section 4.2). 

Also models in which a leptonic global symmetry group G is spontaneously broken 
into a conserved lepton number L, G ^ U~{1), have been considered, and the ensuing 
Goldstone bosons are also called Majorons many times. In this case the Majorons 
carry the conserved L number. L may be a non-orthodox lepton number (such as 
Le + Lfj^ — Lr, for example. The existence of this (almost) conserved lepton number 
produces the appearance of at least one (pseudo) Dirac neutrino mass eigenstate. 
Models of this type were produced for the "17 keV neutrino" and to explain a possible 
signature for neutrinoless double beta decay with emission of a boson (see later). 

Finally, interesting models result when an explicit L-breaking is present in Majoron 
models, transforming the Majoron into a massive pseudo-Golstonc boson. The explicit 
breaking may be due to quantum gravity (Akhmedov, Berezhiani and Senjanovic' 
1992, Akhmedov et al 1993b, Cline et al 1993), although one may protect the global 
L-symmetry with a gauge symmetry (gauge symmetries are respected by gravity, 
Rothstein et al 1993) or by other mechanisms (Lusignoli et al 1990). In this models, 
Majorons could even be the dark matter in the universe (Akhmedov et al 1993a, 
Rothstein et al 1993, Berezinsky and Valle 1993). 

A brief comment on familons is in order. Familons F are the Goldstone bosons 
associated with the spontaneous breaking of an inter-familial global symmetry (Reiss 
1982, Wilczek 1982, Gelmini, Nussinov and Yanagida 1983). The main difference 
between Majorons and familons is that Majorons have usually much larger couplings 
to neutrinos than to other fermions, allowing for faster neutrino decays, for example. 
Since the inter-familial symmetry rotates in principle whole families together, or at 
least whole lepton doublets together (in only leptonic inter-familial symmetries) the 
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same coupling allowing i/^ ^ v^F or i/^- ^ i/gF would allow /i — > eF or r — > eF. From 
experimental bounds in these last modes one gets large lower bounds on the breaking 
scale, of order 10^° GeV and 10^ GeV respectively (D'Ambrosio and Gelmini 1987). 

3.6. Models with Gauged Lepton-Number 

The simplest models of this type arc left-right symmetric models, where right and left 
handed fcrmions are assumed to play identical roles, consequently there is a uji for 
every i/^ (see for example Mohapatra and Pal 1991). The maximal symmetry of these 
models is SUl{2) x SUii{2) xU{l) x P, where P is a sort of parity that exchanges left 
and right indices (P-parity insures also that the coupling constants of the two SU (2) 
groups coincide). In the simplest versions of these models there are only three Higgs 
fields, two triplets Aji = (1, 3, 2) and = (3, 1, 2) and the usual doublet = (2, 2, 0). 
The VEV (A/j) ^ gives masses to the right gauge bosons Zr, and Majorana 
masses to the Ufi of the same order M, breaking the symmetry into the SM. In these 
versions, (A/,) ~ Aw^/ (Ai?) results from the minimization of the potential (Mohapatra 
and Senjanovic 1981) and the mixed mass matrix in (2.12) is 



where mn ~ v and A, / are combinations of coupling constants. Consequently, the 
light neutrino masses are, (see (3.18)), 



and, unless A is extremely small, the first term dominates. These models have been 
used, adding an inter-familial (or "horizontal") symmetry, to obtain three almost 
degenerate mass cigcnstates, whose splittings are due to the see-saw mechanism (in 
models that try to account for most of the present hints for non-zero neutrino masses, 
namely accommodating solar and atmospheric neutrino deficits, hot dark matter and 
giving rise to neutrinoless double beta decay close to the present bounds). 

By complicating the left-right models a bit, P may be broken at a larger scale than 
the rest of the group and (A^) becomes negligible. Then a usual see-saw mechanism 
provides small neutrino masses. Phenomenologically M > 1 TeV, due to the bound on 




(3.26) 



= /(Ai) 




(3.27) 
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additional Z bosons and on Wji bosons. With M ~ 1 TeV one obtains the hierarchy 
of neutrino masses i/g ^ " — MeV. 

The solution to the solar neutrino problem may require much smaller masses, 
(see section 7) 10~^ eVilO""^ eV:l eV, pointing to much larger v^i Majorana masses, 
M ~ 10^° — 10^^ GeV, precisely in the range of the intermediate mass scales in 
50(10) Grand Unified models, necessary to obtain the unification of the SM coupling 
constants after the precision measurements at LEP (Mohapatra and Parida 1992, Babu 
and Mohapatra 1993). These precise measurements of the coupling constants at the 
electroweak energy scale showed that these couplings do not converge simultaneously 
to a single value at a large energy scale, the Grand Unification scale M[/, (as necessary 
in Grand Unified models, Langacker 1981, Mohapatra 1986) in non-supersymmetric 
models with one stage of unification. Thus either supersymmetric models or non- 
supersymmetric models with several stages of partial unification are indicated by these 
measurements (see for example Langacker and Polonsky 1992). 

The natural Grand Unified symmetry that incorporates left-right models is 
S'O(IO). S'O(IO) predicts the existence of the v^i by incorporating all Weyl spinors 
of each generation into a 16-dimensional multiplet. In models with intermediate 
mass scales, SO{VS) breaks into the SM in two stages. The first, at the scale 
Mu ~ 10^^ GeV, brings the symmetry to a left-right symmetric group, that breaks at 
the intermediate scale Mj into the SM. In the intermediate stage there is a discrete 
symmetry, called -D-parity, that has the same role of the P-parity mentioned above. 
D-parity needs to be broken at a scale larger than M/ to obtain a see-saw mechanism 
with M ~ M/, otherwise a solution like (3.27) is obtained (Caldwell and Mohapatra 
1994, loannisian and Valle 1994, Bamert and Burgess 1994, Joshipura 1994, Lee and 
Mohapatra 1994). In 50(10) models, mo are similar to the masses of up-quarks. 

Supersymmetric 50(10) also yields unification of the three coupling constants of 
the SM and in this case a see-saw model with My ~ 10^^ GeV produces masses 
10~^^ eV:10~^ eV:10~^ eV. But also in supersymmetric 50(10) intermediate scales 
can be "gravity induced" (Cvetic and Langacker 1992). There are still many other 
possibilities, using other grand unified groups like E^^ technicolour, etc.. 
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4. Neutrino Mass Searches 



4.1. Direct Searches of Neutrino Mass 

These searches are based on kinematical arguments and assume only finite neutrino 
masses. (For a detailed review on direct searches until 1988, see Robertson and 
Knapp 1988). Indirect searches require additional criteria such as Majorana masses 
in neutrinoless double beta decay and large enough mixing angles in oscillation 
experiments (that measure actually not masses, but mass differences). Neutrinos 
were originally proposed by Pauli in 1932, to account for the continuous spectrum of 
the emitted electrons in weak nuclear /3-decays, {A, Z) — > (A, Z + 1) e~u^. If only e~ 
were emitted, their spectrum would be a 5-function with energy equal to the Q-value, 
Mz — Mz+i, the difference in mass between the initial and final nuclei. The minimum 
energy of is its mass, m,^^ (CPT insures u^l and ivehY have the same mass). Thus, 
the upper end-point Eq of the electron spectrum Eq = Mz+i — Mz — tUi,^, and what 
is more relevant experimentally, the curvature of the spectrum near the end-point, are 
sensitive to m,^^ . The best constraints on rriu^ come from the shape of the spectrum 
of electrons emitted in the Tritium decay, ^He e~i'^, because the low Q-value of 
this decay. Eg = 18.58 keV, allows to notice the effect in the spectrum of relatively 
smaller 171,^^. The first bound using this method, m,y^ < 1 keV, was obtained as early 
as 1948 (Curran et al 1949), and by the 1970's the limit was already 55 eV (Bergkvist 
1972). The claim of a positive detection in 1980 by Lyubimov et al originated many 
new experiments attempting to verify this result (for a review see Holzschuh 1992). 
Six of these experiments (see table 2) have now upper limits that rule out the revised 
Moscow claim of 17 eV < rut,^ < 40 eV (Boris et al 1987). 

A striking feature of these experiments can be seen in the second column of the 
table 2, namely they all find negative niy^ values. This means the experiments find 
the e~ spectra near the end-point deformed with opposite curvature with respect 
to what a neutrino mass would cause. They use then a statistical analysis method 
(prescribed by the Particle Data Group 1986) for dealing with a non-physical result, to 
take into account only positive m^^ and obtain the quoted bounds on rrii,^. However, 
the negative m^^ values are large enough that statistical fiuctuations as a cause have 
a very low probability (about 1%). It seems at this point that the systematics of these 
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experiments are not well understood (Wilkerson 1993). Thus, even if nominally the 
combined results of the first five entries of the table 2 imply m,y^ < 5 eV at 95%C.L., 
due to the systematic uncertainties it is hard to claim more than an upper limit of 
about 10 eV. The still unpublished new result from Troitsk involves a new kind of 
systematics and its low negative seems promising. 



Table 2. Present upper bounds on m,y^. 



Experiment 


ml^ ± CFstat ± o-syst[eV^] 


95 % CL hmit 
on rriy^ [eV] 


Reference 


Los Alamos 


-147 ±68 ±41 


9.3 


Robertson et al 1991 


INS, Tokyo 


-65 ± 85 ± 65 


13 


Kawakami et al 1991 


Zurich 


-24 ± 48 ± 61 


11 


Holzschuh et al 1992 


Mainz 


-39 ± 34 ± 15 


7.2 


Weinheimer et al 1993 


Livermore 


-130 ± 20 ± 15 


8 


Robertson 1994 


Troitsk 


-18 ±6 


4.5 


Belesev et al 1994 



The best bounds on m^^ come from the measurement of the momentum of 
produced in the decay of pions at rest, tt iiv^ (Abela et al 1984, Jeckelmann et 
al 1986). After several recent objections (Wilkerson 1993) have been cleared up, the 
present bound is nii,^ < 160 keV at 90% C. L. (Assamagan et al 1994). 

The Ur has only been observed as missing energy in r-decays, contrary to the 
other two neutrinos that have also been observed by producing their corresponding 
charged lepton in interactions with matter. Rare multi-particle semileptonic decays 
of the T provide the best bounds on Ur (Weinstein and Stroynowski). The ARGUS 
collaboration (Albrech et al 1992) obtained m,^^ < 31 MeV at 95%C.L., with 20 events, 
studying the reaction t~ — > 57r Ur- The CLEO collaboration (Cinabro et al 1993), 
with a much larger sample of 113 events, obtained m^^ < 32.6 MeV at 95%C.L.. There 
is already a preliminary upper bound from Beijing of 29 MeV at 95%C.L. (Darriulat 
1994). Improvements on this method are not expected to lower the niy^ bound to 
better than about 10 MeV. 

Besides these laboratory searches, several cosmological and astrophysical 
arguments considered later in this review (sections 5 and 9) provide additional direct 
bounds on neutrino masses. 
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Here we have called m^^ , rrii,^ and rrti,^ the three mass eigenvalues (more properly- 
called 777-1, "72 and 7773) of the mass eigenstates vi, and 1^3, and we assume vx consists 
mainly of v^^^ of and v-^ of v^. In recent years evidences for a 17 keV neutrino 
seen in nuclear beta decay spectra were claimed by various experiments, giving origin 
to a quite intense experimental and theoretical activity (for a review see, for example, 
Gelmini, Nussinov and Peccei 1992). This neutrino was seen as a shoulder in the e~ 
spectrum, at 17 keV of the end-point, indicating that v^^ was composed mainly of a 
very light (or massless) neutrino mass eigenstate with an admixture of order 1% of 
the heavy neutrino mass eigenstate. The existence of this shoulder has been definitely 
rejected in several recent experiments (Hime 1993, Bonvicini 1993). 

4.2. Double Beta Decay 

Neutrinoless double beta decay (PPOu) is a process that provides a very sensitive 
probe of Majorana neutrino masses. Unlike the ordinary double beta decay with 
emission of two neutrinos {f3(32v)^ that is just an allowed but rare transition at second 
order in the weak interactions {A^ Z) — > [A^ Z + 2) + 2e~ + 2v^^ the /J/JOi^ transition 
(^, Z) — > (A, Z + 2) + 2e~ is a process that violates lepton number by two units and 
hence requires a departure from the SM. 

The simplest scenario in which the Pf30u decay takes place is precisely in the 
presence of Majorana neutrino masses, that convert the (i^eL)^, emitted in association 
with one of the electrons, into a that can thus be absorbed in the second vertex. 
Although the required chirality flip introduces a strong suppression in the amplitude 
(with respect to that of the f3f32v)^ the neutrinoless decay has also a larger available 
phase space. Furthermore, it has the much cleaner experimental signature of producing 
a single peak in the spectrum of the sum of the two electron energies. The f3f32v decay 
leads instead to a continuous spectrum that is hard to identify above the experimental 
radioactive background. These considerations also imply that the chances of observing 
a 1313 signal increase if elements with large Q-values are used, because of the significant 
reduction of background at large energies. Another obvious requirement is that the 
single 13 decay of the isotope be absent or strongly suppressed. 

Since the typical lifetimes of emitters are 10^^-10^^ yr, the first indirect 
experimental evidence for this process came from geochemical searches, studying 
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elements for which the daughter nucleus is a noble gas (and hence naturally absent in 
solid materials). Looking for instance for ^^°Xe and -^^^Xe in Te rich rocks or for ^^Kr 
coming from ^^Se, the double beta decays of these Te and Se isotopes were established. 
These experiments however cannot distinguish between the 2v and channels, with 
the possible exception of the ratio of ^^^Te and ^^*^Te lifetimes (Bernatowicz et ai 
1992), due to the similarity of the corresponding matrix elements and the difference of 
phase space ratios in both types of decays. A somewhat similar radiochemical study 
of a 238^ artificial sample revealed the presence of ^^^Pu atoms produced by decay 
after the original purification (Turkevich et al 1991). 

In the last few years, due to the significant improvements in background reduction, 
detector technology and the use of large amounts of isotopically enriched samples, the 
l3/32iy process was directly observed in the isotopes ®^Se (Elliot et al 1992), ^^Ge 
(Balysh et al 1994, Avignone et al 1991), ^^^Mo (Elliot et al 1991, Dassie et al 1994, 
Ejiri et al 1994), ^^^Cd and ^^*^Nd (for recent experimental reviews see Moe 1993 and 
Morales 1992). The measured lifetimes are also in satisfactory agreement with recent 
theoretical calculations (Wu et al 1991). This has given renewed impetus to this field. 

Many experiments are also devoted to pursue the more interesting search of the 
neutrinoless mode. The non observation at present of any peak at the endpoint of the 
two electron spectrum implies lower bounds on the lifetime of the PPOv decay, T°/'2- 
The Heidelberg-Moscow 90% C.L. constraint T^J^C^Ge) > 1.9 x 10^^ yr (Maier 1994), 
provides the most restrictive constraint on the effective Majorana e-neutrino mass (see 
below), (m^^) < 1.1 eV (the precise bound slightly depends on the theoretical matrix 
clement calculation followed). Other bounds have been obtained for instance from 
T^J^i^^^Xe) > 3.4 X 10^3 yr, which implies {m^J < 2.8 eV (Vuilleumier et al 1993). 
In the near future, the experiments with ^^Ge, and probably also those with ^^^Xe, 
loojVIo and ^^^Cd, are expected to reach a sensitivity down to {m^J ~ 0.1 eV. 

If the electron neutrino is a Majorana mass eigenstate, the effective e-neutrino 
mass (rrii,^) which appears in the transition amplitude is just m^^. If the electron 
neutrino is instead a linear combination of several mass eigenstates, I'e = 5^ t^em^^m; 



where we use self-conjugate Majorana fields fm with positive masses rriiy^ . The matrix 



one has 




(4.1) 



32 



Uem is in general complex. If CP is conserved, the arbitrary phases in the fields can be 
chosen so that the matrix U becomes real, but in so doing, if the Majorana fields are 
taken to be self-conjugate their masses may not be all positive, i.e. m^^ = A^|m,y^|, 
with Am = ±1- The effective mass can be written in the CP conserving case as 



(4.2) 



We recall that rjcp = is the intrinsic CP parity of the Majorana neutrino (see 
section 2, Wolfenstein 1981, Kayser 1984, 1988). 

Equation (4.2) shows that cancellations can occur between different states in {nijj^). 
If this happens, the actual neutrino masses can be all larger than the bound on {rrit,^) 
implied by the non-observation of PPOv. In particular, a Dirac can be viewed as the 
limiting case of two degenerate Majorana states with opposite CP parities (see section 
2), for which the cancellation is complete and no neutrinoless decay takes place. This 
is of course also expected from the absence of lepton number violation in this case. If 
CP is violated in the neutrino sector, there will be unremovable phases in the matrix 
U and hence cancellations can occur even among states with equal CP parities. 

Eq. (4.1) is actually valid if all the Uj^ states appearing in z/g are light {m^^^ < 10 
MeV), since otherwise the nuclear matrix elements involving the heavy neutrino 
propagator are strongly suppressed. Hence, very heavy components have an additional 
suppression besides that of the small factor. 

In addition to the "mass mechanism" just discussed, the PPOu can also take place 
in left-right symmetric theories. In this case, the presence of vu states, right-handed 
charged currents as well as a larger Higgs sector combine to give rise to the neutrinoless 
decays in several possible ways (for a review see e.g. Vergados 1986, Tomoda 1991). 

Other interesting extensions of the SM affecting /3/3 processes are the Majoron 
models (see section 3.5). Here the coupling of the Majoron Goldstone boson (J) to 
neutrinos 

JC = -i^^^l^ 1^{K)rJ + h.c. (4.3) 

produces the required chirality flip in the exchanged neutrino line by the emission 
of the neutral Majoron*. Due to the energy carried away by the undetectable J, 

* Since LEP has excluded both triplet and doublet Majorons, this effective coupling can only be with 
singlet or mixed Majorons and could arise from a fe mixing with singlet neutrinos (see section 3.5). 
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the electron spectrum is no longer a peak. PPOi/J is however still concentrated at 
larger energies than the ^^2^ spectrum. The non-observation of this decay has then 
resulted in significant constraints on the coupling {geii)u^v^j- In particular, the best 
direct bound Qeff < 7 x 10~^ comes from -"^^^Nd (Moe 1994), even if the lifetime limit 
is weaker than that from Ge (Beck et al 1993), due to the favorable matrix element 
and phase space factors of the first. 



5. Neutrinos in Cosmology 

Cosmological and astrophysical arguments are a complement to laboratory 
experiments as a probe of neutrino physics and actually provide some of the most 
restrictive constraints on certain neutrino properties. 

The possible cosmological consequences of neutrinos, and in particular of their 
non-vanishing masses, started to be explored in the seventies. Bounds on the masses 
from the contribution of stable neutrinos to the density of the universe were derived 
first. Then, bounds on unstable neutrinos were obtained from the effects of the decay 
products. These cosmological tests derived for neutrinos are now routinely applied to 
any new proposed particle to limit its lifetime, mass, cosmological density and decay 
modes. 

The hot Bing Bang, the standard model of cosmology, establishes that the universe 
is expanding from a state of extremely high temperature and density. The moment 
the expansion started is taken as the origin of the lifetime of the universe t, which 
is now determined to be to ~ 1.3 to 1.7 x 10"'^° yr (from the oldest globular clusters, 
Demarque, Deliyannis and Sarajedini 1991 give 14-17 Gyr, Renzini 1993 gives 13-15 
±3 Gyr). Notice that in cosmology the subscript zero denotes present values. The 
Big Bang model is based on three major empirical pieces of evidence, namely: the 
Hubble expansion, the cosmic blackbody microwave background radiation, CMB, and 
the relative abundance of the light elements (up to^Li). These pieces of evidence are 
very difficult to explain with a model different than the Big Bang (Peebles et al 1991). 

The Hubble parameter provides the proportionality between the velocity of 
recession of far away objects v, and their relative distance d, v = Hd. Its present 
value, the Hubble constant, is known up to a factor of two, Hq = lOOh km/sec Mpc, 

34 



with 0.4 ^ h ^ 1, where a parsec is pc = 3.26 hght years (see e.g. Jacoby et al 1992; 
Fukugita, Hogan and Peebles 1993; Scott et al 1994). 

The CMB was produced at t^ec — 3 x 10^ yr, the recombination epoch, when atoms 
became stable and replaced ions and electrons in a plasma as the constituents of the 
universe. The Cosmic Background Explorer (COBE) satellite (Mather et al 1994) 
confirmed that the CMB has a blackbody spectrum, and made the most accurate 
measurement of its temperature. To = 2.726 ± 0.010 K (95% C.L.). Thus we know 
with great accuracy the number density of relic photons (the CMB photons are the 
most abundant in the universe, by several orders of magnitude, see e.g. Kolb and 
Turner 1990, p. 143) = {2C{?>) / tt'^)T^ = 411 cm-^. 

The cosmological abundance of ^He and of the trace elements D, "^He and ^Li are 
well accounted for in terms of nuclear reactions that occurred at t ~ 10~^ — 10^ 
sec, the nucleosynthesis epoch, when the photon temperature was T ~ 10 — 0.1 
MeV (necessarily below the binding energy of the light nuclei). This is the earliest 
available proof of the consistency of the Big Bang model. The relative amount of 
the primordial light elements synthesized is very sensitive to the relative abundance 
of baryons (protons and neutrons) over photons, r] = ub/tIj (n is number density), 
and to the expansion rate (and, thus, to the energy density of the universe p, since 
H ~ y/p) at the epoch of nucleosynthesis. Recent analysis (Walker et al 1991, Smith 
et al 1993) find that r) must be in a smaU range rj = (2.8-4) x 10"^^ (95%C.L.) to fix 
three very different abundances: ^He, D + ^He and ^Li. This implies a bound on the 
energy density of baryons, = Pb/Pcj through the relation rj — 2.68 x 10~^QBh^ , 
which is 0.01 < < 0.015. It is convenient to express energy densities p in units of 
the critical density pc = SH^/SnG = 10.5/i^ keV cm~^ (defined to be such that for 
Ptotal < Pc the universe will expand forever and for Ptotal > Pc the universe will 
eventually recoUapse), i.e. 0, = p/pc- Attempts to avoid the nucleosynthesis upper 
bound on Q,B have been largely unsuccessful so far (Malaney and Mathews 1993). 
These include, inhomogeneities in the baryon density (Applegate, Hogan and Scherrer 
1987), late decaying particles (Dimopoulos et al 1988) and generating entropy after 
nucleosynthesis (Bartlett and Hall 1991). We mention later a new idea involving a 
decaying Uj. (Gyuk and Turner 1994). 

An increase in the expansion rate of the universe leads to an earlier freeze out 
of the weak interactions that determine the ratio of neutrons over protons and. 
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consequently, to a larger value of this ratio and to overproduction of ^He (since as 
a first approximation all neutrons end up in ^He nuclei). Thus, the observational 
upper bound on the primordial ^He abundance provides an upper bound on the total 
p at the moment of nucleosynthesis, that is expressed in terms of the allowed number 
of light neutrino families, N^^ ^ 3 + ^A'',^. Walker et al 1991 found dN^, = 0.3, but 
this number changes often, with the incorporation of new measurements of the cosmic 
abundances of elements*. Even if the total width of the Z° measured at LEP insures 
that Ni^ = 3, the nucleosynthesis bound is still useful. This is because it limits any 
extra contribution beyond those of three left handed relativistic neutrinos to p during 
nucleosynthesis, to be equivalent at most to 5Ni, of a neutrino species. 

Both CMB and nucleosynthesis provide a plethora of bounds on neutrinos, some 
of which we will mention below. Other bounds are provided by the present cosmic 
energy density and by structure formation in the universe (namely the formation of 
galaxies and galaxy aggregates). 

May be the most important cosmological constraint on stable neutrinos is the mass 
bound. Light neutrinos (rrii, < 1 MeV) with SM interactions are kept in equilibrium 
with charged leptons and photons in the cosmic plasma (due to weak interactions) 
until a temperature T ~ 1 MeV. At these temperatures the rate of weak processes 
becomes smaller than the expansion rate of the universe and neutrinos "decouple" or 
"freeze out". Their number (per comoving volume, i.e. a volume increasing due to 
the Hubble expansion) becomes constant afterwards. Taking into account that the 
number of photons is increased when e"'"e~ annihilate (at T ^ mg, when e"'"e~ can 
no longer be formed), due to entropy conservation, the present number density of 
neutrinos (plus antineutrinos) per species is rii.^ = (3/ll)n^ = 102 cm~^. 

If neutrinos are massive, their energy density is Pl,^ = mi/^rii,.^, therefore (Gerstein 

* Loopholes of this bound exist but they require rare specific neutrino properties, such as excess 
lepton numbers of the order of the photon number density n~^, namely {riu^ — n^c') jn-^f ~ 0(1), (while 
one expects this leptonic asymmetry to be of the same order of magnitude of the baryonic asymmetry 
f] ~ 0(10^"'^'') (Olive et ai 1991), or a heavy tau neutrino with mass in the MeV range, decaying 
with a short lifetime ( ^ 30 sec) into v^^ or, better, and other invisible particles (a Majoron, for 
example), as mentioned later in this section (Kawasaki et ai 1994 and Dodelson, Gyuk and Turner 
1994a). 
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and Zeldovich 1972, Cowsik and McClelland 1972) 

i=l 

This bound assumes that only left-handed neutrinos are ever in equilibrium in the 
primordial plasma (what is correct, even in the case of Dirac neutrinos for Dirac 
masses below a few keV, see e.g. Kolb et al 1991). The best upper bound on Qoh^ 
docs not come from estimates of Oq and h, but from to- Due to the relation between 
to, Ho and Qo (that depends on the nature of the content of the universe: matter, 
radiation or a cosmological constant, namely vacuum energy), a lower bound on the 
age of the universe translates into an upper bound on Qoh^- Thus, to ^ l-Sx 10^^ yr 
provides the bound flgh'^ ^ 0.4 for a matter dominated universe (and zero cosmological 

3 

constant, see Kolb and Turner 1990, figure 3.3), thus ^ niiy. ^ 37 eV. Only for small 

i=l 

values of h, 0.4 < h < 0.5, it can be 0.25 ^ Qoh'^ ^ 0.4 (see the same figure). 
Thus, if the more popular lower bound of /i ~ 0.5 is taken then JIq/i^ ^ 0.25 and 

3 

consequently, ^ m^. ^ 23 eV. The presence of a cosmological constant A would relax 

i=l 

these constraints a bit (because Cloh'^ can be larger given the same to, see e.g. Kolb 
and Turner 1990). Neutrinos with mass close to these upper bounds could dominate 
the mass density and provide, therefore, the Dark Matter (DM). 

It is by now well established that the dominant form of matter in the universe 
is only detectable through its gravitational effects (and, because of this, it is called 
Dark Matter). Most measurements of the DM at large scales obtain ^dm — 0.2 
to 0.3, but some obtain close to 1 (see e.g. Kolb and Turner 1990, Peebles 1993). 
Because all other known contributions to Qq are much smaller, the total amount of 
DM in the universe is responsible for the ultimate fate of our universe, expansion 
forever or recoUapse (if < 1 or > 1, respectively). The nature of the DM is 
one of the most important entirely open questions in physics. DM candidates are 
classified as cold or hot according to the galaxy formation scenarios derived from 
them. Galaxies are assumed to be formed through gravitational instability, from tiny 
inhomogeneities in the energy density {6p/p < 10~^ at recombination), that leave an 
imprint in the CMB, recently found by COBE. Very different mechanisms result if the 
DM is relativistic (hot DM, HDM) or non-relativistic (cold DM, CDM) when galaxy 
size inhomogeneities could first start collapsing at T ~ 1 keV. Massive neutrinos (with 
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m < 1 keV) could be HDM. Simulations of structure formation with pure HDM fail 
to fit the data, because galaxies form too late (by fragmentation of the much larger 
structures that form first). 

Pure CDM accounts for the bulk of the known data, but seems to fail in detail. The 
COBE measurements of anisotropies in the CMB provide a measurement of density 
inhomogeneities (Smoot et al 1992, Gorski et al 1994). In the context of models with 
fl — 1 and primordial scale-invariant fluctuations (the simplest assumptions), once 
the normalization given by COBE is imposed on the spectrum of density fluctuations 
predicted by pure CDM at large scales, the spectrum has too much power on smaller 
scales, namely the scales of galaxy clusters (Wright et al 1994). A possible acceptable 
solution explored at present is mixed DM (MDM or HCDM). Recent simulations 
suggest an admixture of ~ 0.30 of HDM in neutrinos in a universe dominated by 
CDM (Wright et al 1994, Nolthenius et al 1993, Bonometto et al 1993, Klypin et al 

3 

1994). These simulations use Qo = 1, and h — 0.5, thus requiring J2 '^v, — ^ ^V. 

i=l 

Up to now we mentioned neutrinos lighter than 1 MeV. Neutrinos with ^ 1 
MeV would decouple while they are non-relativistic and their density is thus reduced 
by a Boltzman factor Q.^, ~ exp(— mj^/T). In this case the bound in the energy density 
Q.oh?' is satisfied for rrii, ^ few GeV (Lee and Weinberg 1972, Hut 1977, Sato and 
Kobayashi 1977, Vysotsky, Dolgov and Zeldovich 1977). The mass bounds differ for 
a Dirac or a Majorana neutrino, they are higher for the latter (Kolb and Olive 1986, 
see also Kolb and Turner 1990). Only a A^^ generation neutrino could be that heavy, 
but the LEP bound excludes its existence unless rrii, > Mz/2 ~ 45 GeV. These heavy 
neutrinos would only have a small density Qi, < 10~^. 

These constraints apply to neutrinos with only standard interactions (i.e. to 
neutrinos that are non-standard only because they have a non-zero mass). Neutrino 
masses could be in the forbidden range mentioned above, namely 30 eV < rui, < few 
GeV, if they have interactions beyond the SM that allow for faster annihilation or 
decay. 

Unstable neutrinos with m < 1 MeV whose relativistic decay products dominate 
the energy density of the universe until the present must have a lifetime 

T < (92 eY /m^fiVtoh^fto , (5.2) 

to insure that the energy density of the decay products, JIdp, is not too large, i.e. 
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^DP < (Dicus, Kolb and Teplitz 1977, Pal 1983, Kolb 1986)*. The bound in (5.2) 
and the corresponding bound for heavier masses are shown in figure 1 with a continuous 
contour fine. Since neutrinos and photons are almost in equal numbers, the neutrino 
energy density = nvTrit, becomes dominant over the radiation density prad — n-yT as 
soon as rrii, > T, i.e. neutrinos matter-dominate the energy density of the universe as 
soon as they become non relativistic (actually at T~ OArriy). Thus, if neutrinos decay 
while non- relativistic into relativistic decay products, i.e. p^, = pop t = r^,, these 
products radiation-dominate the universe. Because in (5.2) the universe is assumed 
to be radiation dominated until the present, the constraint to > 1 x 10-*^° yr (1.3 
xlO-'^° yr) requires Q.oh? < 0.3 (0.1). If the universe is assumed to become dominated 
by matter at some time after the neutrino decay, the bound on the lifetime is more 
restrictive than (5.2) (and the bound in figure 1)**. In fact, the massive neutrino 
has to decay earlier for the energy of their decay products to have a longer time to 
decrease, as T^, and become subdominant before the present with respect to a matter 
density component that decreases slower, as . This is actually required by structure 
formation arguments. Because the growth of density fluctuations is suppressed in the 
period of radiation domination of the decay products, this period should finish at 
most when structure in ordinary matter could start forming, i.e. at recombination 
tree ^ 10~^to, when atoms become stable. This argument replaces to by tree in (5.2), 
thus yielding a bound more stringent by a factor 10~^ (Steigman and Turner 1985)***. 
This bound is shown for m < 1 MeV in figure 1 with a dashed contour line. 

A neutrino with mass in the keV to MeV range, with a lifetime close to the lower 
limit just obtained, could not only be harmless but actually can help in the structure 
formation. It could actually be what CDM needs to find perfect agreement with 
data (White, Gelmini and Silk 1994, Bardeen, Bond and Efstathiou 1987, Bond and 
Efstathiou 1991). The main effect of the decaying neutrino is to delay the onset of the 

* This bound is obtained with a simple calculation that assumes all neutrinos decay suddenly at 
t=T. A proper calculation shows that this assumption overestimates the Odp by about 15% (Turner 
1985). Thus the bound (5.2) should be actually higher by a factor of 1.3. 

** Notice this is in contradiction with the equivalent bounds derived in Kolb and Turner 1990 and 
presented in their figure 5.4, that we believe are not correct. 

*** For a more accurate version of this bound based on numerical calculations and using the present 
constraints from large scale structure and CMB see White, Gelmini and Silk 1994. 
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matter domination by the CDM, by adding to the radiation density the contribution 
of the (relativistic) decay products. Also a heavier neutrino, a u^- with mass between 
1 and 10 MeV could do it (Dodelson, Gyuk and Turner 1994b), provided the lifetime 
is ^ 10^ sec to avoid excessively perturbing nucleosynthesis (Dodelson, Gyuk and 
Turner 1994a, Kawasaki et al 1994) 

In fact, a non-relativistic neutrino present during nucleosynthesis, necessarily a 
tau-neutrino, may contribute more to p than a relativistic species. Therefore, the 
bound on the additional number of relativistic neutrinos mentioned above, SNi, < 0.3, 
forbids the mass ranges 0.2 MeV-33 MeV (for a Dirac neutrino) or 0.4 MeV-30 
MeV (for a Majorana neutrino), if the neutrino is present during nucleosynthesis, i.e. 
for T > 10^ sec (Dodelson, Gyuk and Turner 1994a, Kawasaki et al 1994, Dolgov, 
Kainulainen and Rothstein 1994, Kolb and Scherrer 1992, Kolb et al 1991). This 
bound is shown in figure 1 with a dot-dashed contour line. These nucleosynthesis 
bound combined with the laboratory upper limits (31 MeV, see section 4) nearly 
exclude* a Uj- more massive than about 0.4 MeV, if r > 10^ sec (and if the Uj- docs not 
have additional annihilation channels besides those of the SM, so that its relic density 
is the one computed with standard interactions). 

For shorter lifetimes the bounds depend on the decay channel. Radiative decays, 
i.e. into photons or e+e", are excluded for all allowed neutrino masses, as can be seen 
in figure 2. Bounds come from not allowing distortions of the CMB (this excludes the 
gray region in figure 2, taken from Kolb and Turner 1990) and from excessive entropy 
generation or the dissociation by the decay products of the synthesized nuclei after 
nucleosynthesis. Also stringent bounds come from astrophysics (section 9). Bounds 
obtained from the supernova SN1987A exclude the hatched region in figure 2. 

The decay modes available are Ur 3z/' and Uj- — > Ue(p. The mode Ut — > UgCp, 
has been studied for all masses and lifetimes. For very short lifetimes (f) is brought 
into equilibrium by decays and inverse decays and its contribution to p during 
nucleosynthesis becomes unacceptably large (Kawasaki et al 1994). 

Let us mention two curious properties of these massive, fast decaying tau neutrinos. 
A Ut- with mass m >1 MeV and short lifetimes r < 10 sec, that decays into another 
neutrino (and other sterile particles), could decrease the amount of ^He produced thus 

* The preliminary bound from Beijing of 29 MeV, if confirmed, would already have closed the gap 
between both limits. 
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actually loosening the nucleosynthesis bound on new light species. Available analyses 
differ slightly in the range at stake (Dodelson, Gyuk and Turner 1994a, Kawasaki et 
al 1994). This possibility could become particularly useful if the observational data 
would lead to a lower ^He abundance than now assumed (reducing the allowed number 
of effective neutrino families to less than three). The second property refers to a Uj- 
with mass (20-30) MeV decaying into Ug (and other invisible particles) with a lifetime 
of (200-1000) sec, that has a relic density smaller than that of a standard neutrino 
(i.e. this neutrino must have larger than standard annihilation cross sections in the 
early universe). In this case the bound on Qb from nucleosynthesis could be loosen 
(Gyuk and Turner 1994). Moreover there is the possibility we already mentioned of 
a massive short-lived neutrino of m > 1 MeV helping in structure formation in the 
early universe (Dodelson, Gyuk and Turner 1994b). Rejecting experimentally a Vt- 
mass larger than 1 MeV would eliminate these possibilities. 

As already mentioned, radiative decays {u z/'7, u z/'e"'"e~) cannot help in 
evading the cosmological mass bound (see figure 2), i.e. cannot have shorter lifetimes 
than those excluded in figure 1. Only decays into neutral weakly interacting particles 
are allowed, namely u — > Su' and u — > v'cf). It is difficult to find extensions of the SM in 
which the first decay is fast enough while avoiding conflicts with cosmological as well as 
laboratory bounds and suppressing the related forbidden radiative decays (for a review 
see e.g. Mohapatra and Pal 1991). The mode v v'cj), where is a Goldstone boson, 
seems the most promising, and finding a neutrino mass in the forbidden cosmological 
range would strongly suggest the existence of Goldstone bosons. 

If neutrinos are "cosmologically stable", i.e. t to, the radiative decay of a 
massive neutrino is not excluded by the above mentioned bounds. A decaying DM 
neutrino with rrii, ~ 0(10 eV) would however produce a monochromatic UV line (with 

= m,^/2 if the daughter neutrino is massless) that could be observable with satellite 
detectors. A neutrino of m,^ ~ 28 eV and r ~ 10^^ sec has even been proposed to 
account for the observed ionization of the galactic and intergalactic hydrogen (De 
Riijula and Glashow 1980, Sciama 1990). Although this lifetime is much shorter than 
those resulting in the SM with the addition of Dirac neutrino masses, it is surprising 
that it is just at the right value in models in which both masses and radiative decays 
(through a dipole transition) are generated by loop effects (Roulet and Tommasini 
1991, Gabbiani et ai 1991). 
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The nucleosynthesis bound on the effective number of relativistic neutrinos SN^, 
apphes also to any new mechanisms that could bring sterile particles into the 
primordial plasma, such as active-sterile neutrino mixings (see sections 7 and 8) or 
a Dirac neutrino mass. With respect to Dirac masses, recall from section 2 that if 
neutrinos are Dirac particles, there exist right-handed chirality neutrinos (z/^j, and 
also left-handed chirality (z^^)l) that in the SM do not have weak interactions. Since 
the physical states of a neutrino are helicity states, and helicity and chirality do not 
coincide for massive neutrinos, all four helicity states have weak interactions. However, 
the admixture of the "wrong" chirality (the interacting one) in a right-handed helicity 
neutrino is of order m/E (for m < E, E is the neutrino energy), thus their interactions 
are suppressed by a factor (m/E)'^ with respect to those of a left handed helicity 
neutrino. An upper bound on the abundance of right-handed helicity neutrinos (given 
by dN,y in this case) translates into an upper bound on their reaction rates and thus 
on the Dirac mass. If m < 0.3 MeV the right-handed neutrinos decouple before 
the quark-hadron QCD phase transition, in which there is a large increase of the 
number of still interacting particles. Thus the equilibrium number of right-handed 
neutrinos becomes relatively very small (Fuller and Malaney 1991, Enqvist and Uibo 
1993). However, there are also out of equilibrium processes in which right handed 
neutrinos are produced, such as pion decays (Lam and Ng 1991). Dolgov, Kainulainen 
and Rothstein (1994) find m^^ ^ 170 keV [{5N^ - 0.10)/0.20] ^2 and m^^ ^ 210 
keV[(5A^j^ — 0.10)/0.20]"'^/^, if the temperature of the quark-hadron phase transition is 
Tqcd = 100 MeV. These bounds are anyhow less restrictive than the bound derived 
from the supernova SN1987A, that forbids Dirac masses from 0(10 keV) to 1 MeV 
(see section 9). Notice that if the right handed neutrino would be trapped inside the 
supernova by unknown interactions beyond the SM (Babu, Mohapatra and Rothstein 
1991) thus avoiding the SN1987A bound (based on energy loss due to the escape 
of right handed neutrinos), the same interactions would bring them into equilibrium 
during nucleosynthesis (leading to an unacceptable 5Njy=l). This conclusion is difficult 
to obviate (Babu, Mohapatra and Rothstein 1992). 

Let us mention last an upper bound on Majorana neutrino masses due to 
the persistence of a baryon asymmetry generated at temperatures higher than the 
electroweak phase transition, as most baryon generation models assume. The L 
violating processes associated with the Majorana mass in conjunction with B + L 
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violating non-perturbative processes in the SM would erase any existing B asymmetry 
unless the Majorana masses are small enough. A stringent upper bound of 0.1 eV was 
formerly derived, but recent re-considerations of this bound have loosen it considerably 
to around 10 keV (see Cline, Kainulainen and Olive 1993b). 



6. Neutrino oscillations 

6.1. Neutrino mixing 

The phenomenon of neutrino oscillations (Pontecorvo 1958), appears because the 
neutrino current (or interaction or flavour) eigenstates Va, Oi = e,iJ,,T (namely, the 
neutrino states produced in a weak decay in association with a given charged lepton 
flavour, see (3.8)) are generally superpositions of different neutrino mass eigenstates 
Ui, i = 1,2,3, 



Here the unitary matrix U — , where K( is the leptonic analog of the Cabibbo- 
Kobayashi-Maskawa quark matrix Kq (see (3.8) and (3.9)). Due to the dissimilar 
propagation of the neutrino mass eigenstates, the flavour content of the propagating 
neutrino v changes with time. Since the neutrinos are usually detected by means of 
charged current processes, sensitive to the neutrino flavour, this oscillating behaviour 
may be observable. Moreover this is an interference effect, sensitive therefore to very 
small neutrino mass differences. 

6.2. Oscillations in vacuum 

To study the u oscillations, assume that a z/q, flavour eigenstate (namely u^, z/^ 
or Ur) is produced at a; = 0, t = 0. Although the u states are wave-packets with a 
certain spread in momentum, it is sufficient to consider just the plane wave solutions, 
since limitations due to the flnite coherence length of the wave-packets in general 
do not affect the cases of interest (Nussinov 1976, Kayser 1981). We know that the 




(6.1) 
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space-time dependence of a free mass eigenstate lyi^t) of momentum pi and energy 



Vp'i + ™i is 



(6.2) 



It is convenient to take the propagating neutrinos with a common definite momentum 
p (the same conclusion result if we take e.g. a common energy and different momenta 
(Winter 1981) or the wave-packets themselves), so that for the ultra-relativistic 
neutrinos £'j ~ p + m^/{2E) and the neutrino state is then 



uit) = exp[ip{x - t)] 22 UaiG^P ( ) ■ 



(6.3) 



The probability with which the neutrino produced as i^q is converted into after 
travelling a distance x = t results 

2 



=Re ^ UaiU*jU^iU/3jexTp 



2E 



pi 



-X 



(6.4) 



We see that in order to have P{i'a ^p) ^af3 it is necessary that at least two 
neutrinos be non-degenerate. 

These probabilities satisfy some important relations. CPT invariance implies 
P(^'q; I'p) = -P(i^/3 — i^a)j while for two flavour mixing, and only if CP is conserved 
also for mixing with three or more flavours, one has P{uoi vp) = Pii'p I'a) 
(Cabibbo 1978), because these probabilities are related by the replacement U ^ U* iiv 
eq. (6.4). These relations are valid when neutrinos propagate in the vacuum. When 
interactions with matter affect the neutrino propagation they may no longer hold 
because the medium itself is generally not symmetric under CP and CPT. Finally, 
the unitarity of U implies the probability conservation relation 



(6.5) 



We will hereafter consider the case of mixing between just two neutrino flavours, 
and vp, 

( COS sva.9\ I vi 
j y — sin ^ cos 6 J \ 1^2 



(6.6) 
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To be more general, Vfj can be taken here as one of the three known neutrinos or 
eventually as a new light singlet neutrino Vs not participating in the weak interactions 
(sterile species). For a sterile neutrino the oscillation formalism is still valid but the 
predictions for observations by both charged and neutral current processes are affected. 
The probability of neutrino conversion then results 

P{T^a T^f}) = sin^ 29 sin^ —tet^ ' (6-7) 

where Am^ = m| — m\ . Clearly the probability of observing the same flavour in the 
neutrino beam is just P{i'a — i^a) = 1 — P{i^a — 1^/3)- We see that these probabilities 
oscillate with an amplitude proportional to the neutrino mixing factor sin^ 29. Noting 
that 

-^ = 0.8mA^ , (6.8) 

Am2 Am2 [eV^] ^ ' 

we also see that the oscillation length is macroscopic for the typical energies of reactor 
(1 — 10 MeV) or accelerator (10~^ — 10^ GeV) neutrinos, at least for the cosmologically 
allowed values of the masses of stable neutrinos (< 10^ eV). 

When considering an experiment, it is important that the neutrino source is not 
monochromatic but instead has in general a broad energy spectrum. Furthermore, 
both the region of neutrino production and the detector have finite sizes. These facts 
make the oscillatory term in the conversion probability to be usually averaged out, 
leading just to P(z/q i>a) = 1 — sin^2^/2. However, if Am^ is so small that the 
oscillation length becomes much larger than the baseline d between the production 
point and the detector, i.e. Am^ [e^^] <^ E [MeV]/d [m], the oscillations do not have 
enough time to develop and no effect can be seen. Clearly going to larger baselines 
allows to test smaller Am^ values, as long as the reduction in the v flux (oc d~'^) does 
not become the limiting factor. Actually, there is intense activity at present to develop 
long-baseline (10 — 10^ km) experiments using accelerator neutrinos, to increase the 
Am^ sensitivity. 

The main sources of artificial neutrinos used in oscillation experiments are: 

- Ue from P~ decays of fission products in reactors; 

- Fluxes of ly^, and z/g, in comparable amounts and with energies of tens of MeV, 
from decays of stopped tt"*" at low energy accelerators; 

- Beams of v^j, or P^, with a small contamination at the percent level of Ve and Pe; 
from decays in flight of tt and K produced in high energy accelerators. 
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In addition to the experiments using artificial sources, studies of v oscillations 
can be performed using natural sources of neutrinos, such as atmospheric neutrinos 
produced by cosmic rays in the upper atmosphere or solar neutrinos produced in the 
fusion reactions inside the sun. In table 3 we show the typical energies, the distances 
involved and the minimum Am^ values leading to testable oscillations, for the different 
kinds of neutrino sources just mentioned. 



Table 3. Neutrino energies baseline distances d and minimum testable mass- 
square differences Am^ for reactor, accelerator (existing short-baselines and proposed 
long-baselines), atmospheric and solar neutrinos. 





reactor 


accelerator 


atmospheric 


solar 






short-base. long-base. 






E [MeV] 


^ 10 


30 - 10^ 


10^ 


^ 14 


d [m] 


10 - 300 


10^ - 103 10^ - 10^ 


10'^ - 10^ 


10" 


Am^ [eV^] 


10-2 


10-1 10-^ 


10-4 


10-11 



As we will see in sections 7 and 8, studies of solar and atmospheric neutrinos 
provide some indications in favour of neutrino oscillations. On the contrary, the 
searches of oscillations using reactor and accelerator neutrinos have provided up to 
the present no compelling evidence for oscillations, therefore excluding the ranges of 
Am^ and sin^ 2Q explored, as we turn now to discuss. 

6.3. Oscillation experiments 

There are essentially two strategies for detecting neutrino oscillations starting from 
a beam of neutrinos of a given flavour i^a' 

i) to measure the survival probability P(z/a ^a) looking for an eventual reduction 

in the z/q, flux, i.e. the so called disappearance experiments; 

ii) to try to directly observe in the detector the interactions due to a different neutrino 
flavour, in the so called appearance experiments. 

Usually a disappearance experiment has two detectors (or one that can be moved) , 
and the comparison of the measurements nearer to the source with those farther away 
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from it are used to search for any possible reduction in the fa. flux due to oscillations. 
This experimental set-up implies that there is not only a minimum testable value of 
Am^ (corresponding to oscillation lengths much larger than the distance from the 
source to the far detector), but also a maximum testable Am^, since for large enough 
mass differences oscillations would have already been averaged in the nearby detector. 

Appearance experiments look appealing because only a few events above 
background are enough to establish an oscillation pattern, while large statistics are 
required for a meaningful signal in a disappearance experiment. There are however 
situations in which an appearance experiment cannot be performed, like in reactor 
experiments, where the CC interactions of the or r neutrinos are kinematically 
forbidden, or for the study of oscillations into sterile species, that only show up in a 
disappearance experiment. 

6.4. Present situation 

In figure 3 we show the constraints obtained on Am^ and sin^ 29, under the 
assumption of two flavour mixing, arising from the unsuccessful searches of oscillations 
of the type — (figure 3a) and — (figure 36). We have only depicted the 
experiments giving the stronger constraints although several more exist. The bound 
are slightly relaxed if mixing among three fiavours is allowed (Bliimer and Kleinknecht 
1985). 

For oscillations of e-neutrinos, the best bounds on Am^ result from reactor 
experiments because of their small energies. The are here detected by the inverse /? 
reaction PgP e'^n. Since these are disappearance experiments, the resulting bounds 
apply to oscillations into any type of neutrinos, and they actually represent the best 
experimental constraints on Ve-J^r oscillations. At the Goesgen reactor in Switzerland 
(Zacek et al 1986) the fliixes measured at three distances (37.9, 45.9 and 64.7 m) were 
compared among themselves and in addition they were compared to the flux expected 
from the knowledge of the u source and of the detector efficiency (this explains why the 
excluded region extends also to large Am^ values) . The experiment of the Kurchatov 
group (Vidyakin et al 1990), in the site of three nuclear reactors in Moscow, used one 
detector at 57 m from the flrst and second reactors and 231 m from the third. These 
comparatively large distances allowed them to exclude down to Am^ ~ 8 x 10~^ (for 
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maximum mixing) from the comparison of the rates measured with all the different 
possible combinations of reactors on and off. 

Although accelerator experiments have not reached such low values of Am^, 

^ i^e appearance experiments have achieved sensitivities to sin^2^ values much 
smaller (~ 3 x 10""^) than those reached at reactors ( ^ 10~^). The best bounds 
on sin^ 29 are from the BNL experiments E776 (Borodovsky et al 1992) and E734 
(Ahrluder et al 1985) and from the Serpukhov experiment SKAT (Ammosov et al 
1988). Also shown in figure 3a are the constraints from the Los Alamos Meson 
Physics Facility (LAMPF) (Durkin et al 1988) and the bubble chamber BEBC at 
CERN (Angelini et al 1986). 

Regarding the — Vj- oscillations, searches looking for disappearance were 
performed at Fermilab CCFR experiment and by the CDHS (Dydak et al 1984) 
and CHARM (Bergsma et al 1984) collaborations at CERN. The best constraints 
on sin^ 2^^,- were obtained by the Fermilab E531 experiment (Ushida et al 1986) with 
an emulsion detector sensitive to the appearance of a r track. Recently, also the 
CHARM II collaboration (Gruwe et al 1993) was able to constrain the Uj- appearance 
by means of the study of the r — > TrVj- decay mode in their fine-grained calorimeter. 
These results are shown in figure 3b. 

There is also an important astrophysical constraint on Ue — ^ix,t oscillations, to be 
discussed in section 9, arising from the non-disruption of r-processes in supernovae, 
that excludes Am^ values bigger than ~ 4 eV^ for sin^ 29 ^ 10~^. 

6.5. Future prospects 

While one may say that in the past oscillation experiments have been performed 
just looking blindly into the attainable ranges of Am^ and sin^ 29, the situation looks 
quite different for the future. In fact, there are two main physics issues that motivate 
the forthcoming experiments. 

First, there is the dark matter problem, that may be accounted for, at least 
partially, if the heavier of the three neutrinos has a mass in the range between a few 
eV and a few tens of eV (see section 5). This is most likely a Ur, and an analogy with 
the mixing in the quark sector would suggest that its most significant mixing is with 
Vfj, (the nearest generation). Thus, this argument has strongly encouraged — > Vr 
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appearance searches sensitive to tiny mixing angles for Am^ ^ eV^ (Harari 1989). 
Two such experiments, CHORUS and NOMAD (di Leila 1993), are now running at 
CERN. They are expected to take data in the period 1994-1995, and to be sensitive 
down to sin 29 ~ 3 X 10"^. CHORUS (de Jong et aJ 1993) has an emulsion target 
with scintillating fibers to assist in the location of the r tracks while NOMAD selects 
kinematically the r decays based on the correlation between the resulting hadronic 
shower and the missing pt carried by the produced fr- An approved proposal (P803) 
related to the main injector project at Fermilab, using 1 ton of emulsion target, will 
extend the sensitivity to Vr appearance down to sin^ 2^ ~ 6 x 10~^ by the end of the 
decade (Kodama 1990). 

The second motivation comes from the observations performed with atmospheric 
neutrinos that, as will be discussed in the section 8, may suggest a mixing between 

and i>T- (or i>e) with sin^ 26 ^ 0.4 and Am^ somewhere in between 10~^ and 10"'^ 
eV^. The need to clarify this issue has encouraged the realization of long-baseline 
experiments, where a very intense v beam is pointed towards a large far away detector 
(Pantaleone 1990, Bernstein and Parke 1991). Several proposals of this kind have been 
studied, including sending a beam from: Fermilab main injector to SOUDAN (baseline 
710 km) (Allison et al 1991); BNL to three new Cerenkov detectors at 1, 3 and 20 km 
(Mann and Murtagh 1993); CERN SPS to Gran Sasso ICARUS detector (730 km) 
(Revol 1993) or to NESTOR (1680 km) (Resvanis 1993); CERN to Superkamiokande 
(8750 km) (Rubbia 1993); KEK to Superkamiokande (250 km) (Nishikawa 1992). 

Since some of these experiments have detectors with the ability to measure 
separately neutral and charged current events, a strategy somehow in between the 
usual appearance and disappearance technics can be used to study oscillations among 
active neutrino species. This strategy consists in using the NC events as a measure 
of the total flux (including the flavours resulting from the oscillations) and the CC 
events to measure the original flavour content. Another similar approach is to use the 
through-going muons produced in the rock outside the detector as a measure of the 
content of the beam and the events contained inside the detector (NC + CC) as an 
indication of the total v flux. 

If the atmospheric v problem were due to oscillations in the Vg. — channel (a 
possibility, however, that seems now excluded, see section 7), these may also be studied 
using reactor and large detectors at distances ^ 1 km. There is a proposal to use 
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the San Onofre reactor in California with a 12 ton detector at a distance of 1 km from 
the reactor (phase I) and a kton detector at a distance of 10-15 km (phase II) that 
will cover all the relevant range. Another proposal is based on the use of the Perry- 
reactor and a kton detector at the 1MB site that is 13 km apart. Finally, there is a 
proposal to use the Chooz reactor in France with a detector 1 km away. 

Besides the proposed detectors just discussed, there are some ongoing experiments 
that will provide results in the near future. The third phase of the oscillation 
experiments at the Bugey reactor in France may slightly improve the sensitivity of 
the Goesgen experiment. The KARMEN detector at the spallation source ISIS in 
the Rutherford Appleton Laboratory (Drexlin et al 1990) is studying — > Ue and 
J^fj, appearance. A very important improvement in the — > v^. oscillations 

will be achieved with the Liquid Scintillator Neutrino Detector (LSND) at LAMPF 
(Whitehouse et ai 1991). 

The approximate sensitivities expected in the experiments discussed in this section 
arc shown in figures 4a (z^e ^y) and 46 {v^ z/,-), together with the region excluded 
at present (continuous lines). Also indicated are the regions relevant for the dark 
matter problem and for oscillations of atmospheric neutrinos. 



7. Solar Neutrinos and Oscillations in Matter 

7.1. The Solar Neutrino Problem 

The most important natural source of low-energy neutrinos {E^ < 18 MeV) reaching 
the earth, except at the moment of a galactic supernova explosion, is our own sun. 
In fact, 98% of the solar energy production arises from the fusion chain that converts 
four protons into a ^He nucleus, two positrons (by charge conservation) and two (by 
lepton number conservation), releasing about 27 MeV and with a maximum i/ energy 
of 0.42 MeV (for these so called pp neutrinos) *. The proton- proton chain also gives 
rise to the production of the heavier nuclei Li, Be and B, although with negligible 

* Since the solar energy flux at the earth is ~ 8.5 x lO^-"^ MeV cm~^ s^^, the integrated pp neutrino 
flux must be ~ 6 X 10^° cm~^ s~^. 
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effects on the overall energy production. The decays give rise to an important 
continuum neutrino spectrum extending up to ~ 14 MeV, while the ^Be produces, 
by electron capture, two neutrino lines at energies 0.86 MeV (90%) and 0.38 MeV 
(10%). Another line at 1.44 MeV arises homp + e+p d + v (pep) while the ^He+p 
(hep) reaction gives a small continuous spectrum up to 18.7 MeV. Finally, in the 
CNO (carbon, nitrogen and oxygen) chain, essentially a ^^C catalyzes the fusion of 4 
protons into ^He to produce the remaining 2% of the solar energy output. The decays 
of the ^^O, ^^N and ^^F produced in the intermediate steps lead to other continuous 
neutrino spectra. All this is shown in figure 5 (Bahcall and Ulrich 1988). 

Although the pp neutrino flux is well established theoretically, the other ones can 
depend significantly on the solar modelling and the input data used in it. For instance, 
the ^B and ^Be neutrino fluxes are sensitive functions of the solar core temperature 
Tc 

$(«B) oc T^^ , $(^Be) cx T| (7.1) 

(while ^{pp) oc T~^-^, actually decreases with increasing Tc due to the conservation 
of the total solar luminosity). Hence, a change on the radiative opacities (namely the 
cross-sections per gram) that determine the energy transport in the sun, can affect 
these fluxes significantly. The opacities themselves depend on the abundances of heavy 
elements and in particular on that of Fe, that contributes sizably to the number of 
free electrons in the plasma. Making different estimates of the latter affect the ®B 
neutrino predictions by as much as 10%*. Another source of uncertainties comes from 
the nuclear cross sections involved. In particular, the ^Be(p, 7)®B cross section is 
known at the energies that are relevant in the sun from an extrapolation of the values 
measured in the laboratory at higher energies, and this procedure may also cause 
discrepancies among different ^B neutrino fiux estimates (see e.g. Dar and Shaviv 
1994). 

Solar neutrinos have been observed by four experiments up to now, with the 
thresholds indicated in figure 5. The pioneering CI radiochemical experiment (Davis 
et al 1990), based on the reaction ^^Cl + Ar + e {El^'' = 0.81 MeV), has 

detected solar neutrinos for more than 20 years. It is sensitive to the ^B as well as the 

* Another possible way of modifying the predictions of Tc is to consider non-standard solar models, 
e.g. invoking large magnetic fields in the solar core, significant turbulence, etc. 
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^Be, CNO and pep neutrinos. The Kamiokande water Cerenkov detector (Hirata et 
al 1991a-b) has been able to observe in real time the ve scattering, having an analysis 
threshold of 7 MeV that makes this device sensitive only to the ®B neutrinos. The 
directional information of the events makes this detector the first 'neutrino telescope' 
able to obtain a 'picture' of an extraterrestrial neutrino source (the sun). Finally, 
there are two radiochemical gallium experiments, SAGE at Baksan (Abazov et al 
1991) and GALLEX at Gran Sasso (Anselmann et al 1994). They use the reaction 
^^Ga + i^e — ^•^^ Ge + e, whose low threshold (£^^'^''=0.23 MeV) makes them sensitive 
to the pp neutrinos in addition to the more energetic ^Be, ^B and CNO ones. In table 
4 we show the rates measured in these four experiments, as well as the theoretical 
predictions based on the standard solar models of Bahcall and Pinsonneault (1992) 
and Turck-Chieze and Lopes (1993). 

Table 4. Solar neutrino measurements and theoretical expectations within the 
Standard Solar Model of Bahcall and Pinsonneault (1992), SSM-BP, and Turck-Chieze 
and Lopes (1993) SSM-TCL. 



Experiment 


Measurement 


SSM-BP 


SSM-TCL 


37C1 [SNU] 


2.23 ±0.23^ 


8± 1 


6.4 ± 1.4 


Kamioka (i^) 


0.50 ± 0.04 ± 0.06 b 


1±0.14 


0.77 ±0.17 


GALLEX [SNU] 


79 ± 10 ± 6 ^ 


131.51^ 


122.5 ±7 


SAGE [SNU] 


7q+18+5 d 
'^^-16-7 







^ Davis 1992, ^ Suzuki 1993, ° Anselmann et al 1994, Abdurashitov et al 1994. 

The most intriguing feature of these four measurements is that they all give values 
significantly below the theoretical predictions. This constitutes the so called solar 
neutrino problem. The explanation of this deficit calls either for modifications in 
the solar models or for new properties of the neutrinos, such as mixings or magnetic 
moments, that could convert the electron neutrinos into less detectable species in their 
journey from the center of the sun to the earth (radiochemical experiments are only 
sensitive to Ug, while the f — e cross section relevant for Cerenkov detection is 7 times 
smaller for i/f^^T- than for Vg, and vanishes for sterile species). Also important is that 
the observed Ga rates are not less than the ~ 70 SNU expected (reliably) from the pp 
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neutrinos alone, since this would otherwise exclude any explanation not invoking new 
neutrino properties. However, if all experimental results are taken at face value, one 
should note that an explanation based on a reduction of the core temperature is not 
possible, since it would reduce much more the flux than the ^Be one, contrarily to 
what results from the comparison of Kamiokande with CI data. 

Turning to the explanation in terms of new neutrino properties, the simplest one is 
to invoke oscillations in vacuum of the Vf, (Gribov and Pontecorvo 1967). In view of the 
typical V energies (~ MeV) and the earth-sun distance (1 A.U. ~ 1.5 x lO"*^^ m), values 
of Am^ ^ lO"-'^-'^ eV^ are required for a significant v^. suppression (see eq. (6.10) and 
(6.11)). However, for Am^ ^ 10~^° eV^ oscillations are completely averaged. In this 
case, the reduction factor results the same in all experiments, and it is larger than 1 /2 
for the mixing among two flavours. It is only for 10~^^ eV^ < Am^ < 10~^° eV^ that 
the oscillation length is of the order of 1 A.U. so that neutrinos of different energies are 
converted in different amounts. With this 'just so' oscillations (Barger, Phillips and 
Whisnant 1981, Glashow and Krauss 1987), that also require sizeable mixing angles, 
it becomes then possible to fit simultaneously all three types of experiments (Barger 
et al 1992, Krastev and Petcov 1992). 

A more beautiful possibility is that of explaining the solar neutrino problem using 
the enhancement of the oscillations induced by the effects of the solar medium, that 
we now turn to discuss. 

7.2. Neutrino Oscillations in Matter 

When neutrinos propagate through matter, a crucial differentiation among the 
neutrino flavours appears due to the fact that (unlike or 1/^-), have charged 
current interactions due to W exchange with the electrons. As noted by Wolfenstein 
(1978), these interactions can affect the pattern of neutrino oscillations and the effect 
can be significantly enhanced (Mikheyev and Smirnov 1985) when a resonance crossing 
takes place, in the so called MSW effect (Bethe 1986, Rosen and Gelb 1986, Parke 
1986, Haxton 1986). 

The charged current interaction can be written as Hcc = {^Gf/V2)J^'^J^J (see 
(2.3) and (3.8)) for energies much lower than the W boson mass. Because there are 
electrons in normal matter but not muons or taus, we are interested only in the e 
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i/g interactions. After a Fierz rearrangement to write them as a product of charged 
conserving currents, we get 



G jp 

Hoc = ~ 75)ej^e7M(l - l^)^e 



(7.2) 



For an unpolarized medium at rest the only non-vanishing component of the electron 
current is the temporal part of the vector current, that is nothing but the electron 
density A^e- Hence, an 'effective' potential energy {Hcc) = {ev\Hcc\^^) — V^GpNe 
will affect the propagation phase of the u^,*. Although the neutral currents are 
also present, they can be omitted, since they affect in equal amounts all active 
flavours giving no net effect to the neutrino oscillations (they are relevant however 
for oscillations into sterile species). 

For definiteness we will consider the case of mixing between two active neutrinos, 
i^e and Ua {o( — /J' or r) 




with Rq 




(7.3) 



{s9 = sinO, etc.). The evolution of the flavour content along the neutrino path can be 
simply obtained from the equation 



dx 





(7.4) 



The matrix is 



= - 




{Hcc) 
-{Hcci 



(7.5) 



The first term in the r.h.s. is the usual one, already appearing in vacuum oscillations, 
and the second term arises from the Ue — e coherent forward scattering. We have 
subtracted a piece proportional to the identity matrix (what changes just a common 
phase) to put in a more symmetric form. 

It proves very convenient to define the matter eigenstates as 




— Rf} 




(7.6) 



This can also be thought of as inducing an index of refraction for Vg- 
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where Rq^ diagonalizes the matrix M^, i.e. 



BLM^B,„=1( . |. (7.7) 




Here = Aw? {a - c2eY + s'^2e with a = 2E{Hcc) / ^rr? ■ The matter mixing 
angle 9m entering in Rq^ satisfies 

Even if the vacuum mixing is very small, there will then be maximum mixing in matter 
{9m = 7r/4) in the 'resonance region' corresponding to an electron density such that 
a — c29. In matter with varying density, the width of the resonance corresponds to 
the densities for which \a — c29\ = \s29\. It is important that if we were considering 
antineutrinos, the sign of {Hqc) would be reversed and no resonance would appear 
(for Am^ > 0). Hence, the presence of electrons in the medium (and no positrons) 
results in very different oscillation patterns for neutrinos and antineutrinos. 
The Vm evolution is determined by 



-Aml^E -ij-9. 



i±9m ^m/^E 



with 



d9m. 1 s29 da 



(7.9) 



(7.10) 



dx 2 {a - c29f + s'^29 dx 
We see that the matter states are the propagation eigenstates in a medium with 
constant density (d^^n/da; ~ 0), since no transitions between them can occur in this 
case. When the density of the medium varies along the neutrino path, transitions 
between matter mass eigenstates are induced by a nonzero d9m/dx. This is the case 
in the sun, where the densities fall nearly exponentially in the radial direction. The 
transitions between matter mass eigenstates are usually negligible unless the neutrinos 
are near the resonance layer, for which the diagonal elements in (7.9) are minimum 
and d9m/dx is enhanced. If P = K^^^l^^m)!^ the probability of z/^ conversion 

in the resonance crossing, the averaged probability to detect an electron neutrino that 
has crossed a resonance is 

P.^.^ = l + (^l-P^c29xc29m , (7.11) 
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where 6m is the matter mixing angle corresponding to the point where the was 
produced, while 6 is the vacuum angle*. 

Under the assumption that the electron density varies linearly in the resonance 
layer, the probability of level crossing at resonance is found to be 

P = e-'^^/^ (7.12) 

where the adiabaticity parameter 7 is 

_ ^ Am^ 8^26> 1 

^ ~ 4:E\dem/dx\ 2E c2e \d\nNe/dx\r ^ ' ^ 

(the sub- index r stands for the resonance value). 

In the adiabatic case, i.e. 7^1, the off-diagonal terms in (7.9) can be neglected 
even at resonance, so that P = and the probability of having a u^, after adiabatic 
conversion is P^^u^ = (1 + c29 x c20^)/2. What happens is that if the production 
point is at densities larger than the resonance value, a Ue is mainly the eigenstate 
{9m — 't/2) and remains so during the adiabatic propagation. Thus, it becomes 
essentially when it comes out from the sun. An almost complete flavour conversion 
results. The resonance condition a = c29 can be written as 

Am^c29 ~ — P -^10-^eV^ , (7.14) 

lOMeVpo 

where po — 150 gr/cm^ is the central density of the sun and we have used A^e = YeP/'^p 
(with Yf, = Ne/{Nn + Np) ~ 0.7 — 0.8 in the sun). We then see that the solar neutrinos 
will meet a resonance only for Am^ ^ 10~^ eV^. 

It could be that in the resonance crossing the propagation becomes non-adiabatic 
(namely, 7^1), making the flavour conversion less efficient. This happens when the 
vacuum mixing angle is small so that the resonance becomes very narrow. In the 
extreme non-adiabatic case (7 ^ 1), corresponding to* 

1 dN, 



Am^s'^29 < 2E 



dr 



~ 6 X 10"^ ( — - — 1 eV^ , (7.15) 
* 10 MeV y ' ^ ' 



* This is obtained by projecting at the production point the Ve into matter eigenstates, following 
then the outward propagation in the sun including an eventual jump between matter eigenstates at 
the resonance crossing, and using the fact that the flavour content of those states changes smoothly 
with electron density until they match the vacuum mass eigenstates when they leave the sun. 

* For r > 0.1i?o, Ne oc exp(-lO.5r/i?0). 
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the V flavour content is essentiaUy unchanged in the resonance crossing and no 
enhancement in Piy^. z/q,) is found. 

The survival probabihties depend strongly on the neutrino energy, as is shown in 
figure 6, in which the effect of the conditions of crossing a resonance and that the 
transition not be extremely non-adiabatic are clearly seen. For instance, for s2Q = 0.1 
as in the figure, only the high energy neutrinos [E ^ 10 MeV) would be suppressed 
if Am^ ~ 10""^ eV^ (adiabatic case) while instead for Am^ ~ 10~^ eV^ it would be 
mainly the neutrinos and not the high energy ones to be converted (non-adiabatic 
case). 

An agreement with the deficit observed at the Davis experiment, with adiabatic 
neutrino evolution, is obtained in two regions of the Am^ — s^lQ plane. The first is for 
Am^ ~ 10~^ eV^, if the resonant layer is not too narrow, i.e. slQ ^ 10~^. The second 
is the so-called "large angle solution", s^2Q ^ 0.7 and 10""^ eV^ > Am^ > 10"^ 
eV^, for which the MSW oscillations start approaching the large-mixing vacuum 
oscillations. A non-adiabatic solution exists for Am^s^2^ ~ 10~^'^ eV^, closing a 
'triangle' in the Am^ — s^2d plane. Due to the different energy dependence of the 
reduction factor in the three regimes, experiments with other thresholds, or capable 
of measuring the neutrino spectra, can distinguish among the three solutions. The 
results of the Kamiokande detector in fact disfavor the Am^ ~ 10~^ eV^ adiabatic 
solution, that would deplete too much the high energy neutrinos, while the gallium 
experiments disfavor Am^ values below 10~^ eV^, that would result into an excessive 
conversion of neutrinos. The combination of all the experiments, then, points out 
to two possible solutions, a non-adiabatic one with 8^29 ~ 8 x 10~^ or a large angle 
solution with 5^26* ~ 0.7, both with Am^ ~ 10~^ eV^, as is shown in figure 7, taken 
from Hata and Langacker (1993b) (see also Hata and Langacker 1994). 

For i^e oscillations into sterile species v^^ the solution is slightly different because 
the Vs do not have NC interactions. Hence, the matter effects in the sun as well 
as the cross section for detection in Kamiokande are affected. As a result (Hata 
and Langacker 1994), the large angle solution (that would anyhow also confiict with 
nucleosynthesis (Shi et al 1993)) is absent for — Vg oscillations and the small angle 
solution is slightly shifted. 

There are two interesting cases in which the neutrino propagation through the 
earth may be affected by matter effects. In these cases no MSW effect, namely the 
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crossing of a resonance layer, really takes place since the density varies only slightly, 
from 3 to 5.5 gr/cm^ in the mantle and from 10 to 13 gr/cm^ in the core. Anyhow, 
oscillations of neutrinos that are resonant can be significantly enhanced since the 
matter mixing becomes maximal. The first is the case of the accelerator neutrinos 
in the planned long-baseline experiments and of the atmospheric neutrinos that pass 
through the earth before reaching the detectors. They would be resonant in the mantle 
if Am^ ~ 3 X 10"'^ eV^ for energies ~ 10 GeV. Hence, matter effects can extend the 
sensitivity to smaller s^2^ea values in this mass range (Carlson 1986, Akhmedov, 
Lipari and Lusignoli 1993, Fiorentini and Ricci 1993). However, very long-baselines 
are required to see any effect because the wavelength of the oscillations in matter, 
= 47r£?/A^, becomes 

3 X lO^m 
- tg2^(p/5 gr cm-3) 

in the resonance (taking Yg = 0.5). This is of the order of or larger than the radius of 
the earth. The second case in which matter effects in the earth arc relevant is that of 
the neutrinos observed by Kamiokande {E ~ 10 MeV). The resonance in the earth 
in this case is achieved for Am^ ~ 3 x 10~^ eV^, and can give rise to a day-night effect 
in this mass range if s26 is not too small (if Am is not too large, Baltz and Weneser 
1987, Bouchez et al 1986). The non-observation of this effect (Hirata et al 1991a) 
leads to the exclusion of the region shown in figure 7 (see also Hata and Langacker 
1993a and references therein). 



7.3. Other Solutions 

Many alternative solutions have been proposed to explain the solar neutrino deficit. 
One of them is based on the interaction of the neutrinos with the magnetic fields of 
the sun, that, in the convective, zone increase in strength in periods of high activity 
(Cisneros 1971, Voloshin, Vysotsky and Okun 1986). The resulting effect would be to 
flip the neutrino chirality to produce either a sterile state (with Dirac type moments) 
or just a different active flavour (with Majorana type |AL| = 2 transition moments). 
This solution became especially attractive to account for an anticorrelation with the 
solar activity of the rates, observed in the Davis experiment but not conflrmed by 
Kamiokande. 
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The required magnetic moments are ~ 10~^^ij,b {i^b = e/2me is the Bohr 
magneton). Although these values are not in conflict with direct experimental bounds, 
< 1.08 X 10~^/xb (Review of Particle Properties 1994), they would result in excessive 
stellar cooling (Raffelt 1990), or conflict with SN1987A in the case of Dirac type 
transitions (see e.g. Barbieri and Mohapatra 1988), as we will see in section 9. 

Another difficulty is that it is not easy to flnd models providing magnetic 
moments of the required size keeping the neutrino masses small, since a loop diagram 
contributing to fx gives a mass term when the photon is removed, unless some particular 
symmetries are invoked (Voloshin 1988). Furthermore, transition type spin precession 
is quenched when the mass difference between states is Am^ ^ 10~^ eV^. However, 
matter effects could compensate the mass terms and lead to an enhanced 'spin-flavour' 
precession in a resonance crossing (Akhmedov 1988, Lim and Marciano 1988). Recent 
analysis of the experimental data in this context can be found in Krastev (1993), 
Akhmedov, Lanza and Petcov (1993) and Pulido (1993). 

The usual MSW picture is also modifled in the presence of other non-standard 
properties in the neutrino sector. For instance, if mixing with heavy singlet neutrino 
species occurs, the Z-mediated interactions can become non-universal in flavour. 
Then, it may even become possible to have a resonance for massless neutrinos, by 
compensating the effect of the charged-current with the neutral-current i^e-interactions 
(Valle 1987). However, due to experimental constraints on the weak boson couplings, 
this may only be achievable in media such as neutron stars, but not in the sun. 

More interesting for the solar neutrino problem are new ffavour changing neutrino 
interactions that appear in supersymmetric or GUT extensions of the standard model. 
They may induce neutrino oscillations even for vanishingly small vacuum mixings 
(Roulet 1991), and their resonant ampliflcation can result in sizeable effects, even 
for small (experimentally allowed) couplings. In this type of scenarios, the resonance 
conversion in the massless case may take place in the presence of large new diagonal 
couplings of Ur (Guzzo et al 1991) (for a general discussion see Barger et al (1991)). 
Finally, the possibility that non-universal gravitational interactions (or new flavour- 
dependent long range forces) affect the solar neutrino oscillations (Halprin and Leung 
1991), and the possibility of matter induced neutrino decay into another neutrino and 
a Majoron in the sun (see e.g. Berezhiani and Rossi 1993) have also been considered. 
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7.4. Future Prospects 



Several new solar neutrino detectors are in construction or under study at present. 
The Sudbury Neutrino Observatory in Canada (Aardsma et al 1987), using 1 kton of 
heavy water, will be able to detect, besides the scattering off electrons, the reactions 
Vxd — > priVx and Vgd — > ppe. By observing the n capture, the deuterium disintegration 
{El^ = 2.2 MeV) will allow to measure the NC rates due to all active neutrino flavours. 
The inverse /? decay will be useful to reconstruct the neutrino spectrum accurately 
above ~ 5 MeV. The comparison of the NC and CC rates can give a clear signal 
of oscillations, independently of any solar model input. The spectral shape can help 
to distinguish among diff'erent solutions. This experiment (starting in 1995) and the 
Superkamiokandc water Cercnkov (Totsuka 1990), with 22 kton fiducial mass (starting 
in 1996), will be able to study with large statistics the neutrino fluxes. 

The study of the Be line is the main objective of the Borexino experiment 
(Arpesella et al 1992). It consists of 100 ton fiducial mass of ultrapure liquid scintillator 
with extremely low radioactive background, in order to be able to detect the ue 
scattering in the window 0.25 MeV < E < 0.8 MeV (90% of the signal coming from 
^Be). The events of larger energies will also permit to study the ^B neutrino fluxes. A 
signal that can be useful is the annual modulation of the Be flux due to the eccentricity 
of the earth orbit. For 'just so' oscillations the modulation will differ from the simple 
behaviour. A counting test facility with 2 ton flducial mass is now being built at 
Gran Sasso. Another experiment, Icarus, a liquid argon TPC to be installed at Gran 
Sasso (3 modules of 5 kton each), could allow to study the ^B neutrinos measuring 
separately CC and NC events. A 3 ton prototype built at CERN has performed 
satisfactorily (Revol 1993). Two experiments to detect pp neutrinos are under study 
at present, both helium-based. Heron (Lanou, Maris and Seidel 1987) would detect the 
rotons produced by neutrino interactions in superfluid helium and Hellaz (Arzarello 
et al 1994) would consist of a helium gas TPC at high pressure. 

These experiments, together with increased statistics in the gallium experiments, 
should flrmly establish the nature of the solution to the solar neutrino problem, within 
the next ten years. 
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8. Atmospheric neutrinos 

The cosmic ray nuclei that impinge on the atmosphere are stopped in the upper 
layers by collisions with the air nuclei. These interactions can produce pions and kaons 
that in turn decay giving rise to neutrinos (e.g. tt"*" —>■ /I'^+iy^ and similarly for n~ and 
K^). Also the muons so produced decay after loosing some energy and produce more 
neutrinos via — > + e + (muons with energies below a few GeV decay before 
reaching the ground). The resulting neutrino fluxes on earth are quite significant, 
and were observed many years ago (Krishnaswamy 1971, Reines 1971) by looking at 
horizontal and upgoing muons produced in the rock surrounding deep underground 
detectors (downgoing muons are dominated by those directly produced high in the 
atmosphere by very energetic meson decays). In recent years, it became also possible 
to observe contained events {Ei, < 2 GeV), produced directly by neutrino interactions 
inside large volume detectors (1MB, Kamiokande, Frejus, NUSEX, SOUDAN). 

Atmospheric neutrinos are now an important subject for all underground neutrino 
detectors. They also constitute a background for proton decay searches, and due to 
the long 'baseline' distances involved (10 — 10^ km), they are of particular interest for 
v oscillation studies. 

Several theoretical computations of the resulting u fluxes have been performed, 
both for the low energy neutrinos that lead to contained events (Barr, Gaisser and 
Stanev 1989, Lee and Koh 1990, Honda et ai 1990, Kawasaki and Mizuta 1991, 
Bugaev and Naumov 1989) and for those of higher energies that produce throughgoing 
muons (Volkova 1980, Mitsui, Minorikawa and Komori 1986, Butkevich, Dedenko and 
Zheleznykh 1989). The different computations typically agree on the absolute values of 
the Ve + i^e and V/j, + fluxes only within 20-30%. The reason is that they are affected 
by the poor knowledge of the primary cosmic ray spectrum and composition, the 
uncertainties in the meson production cross section, the use of different calculational 
methods, etc. Also the resulting absolute values of muon fluxes, that could in principle 
be used to normalize the calculations, are poorly known experimentally. Hence, the 
measurements of the absolute neutrino fliixes can hardly give reliable information on 
neutrino oscillations. 

This difficulty can be overcome by considering the ratio R = + v^jive. + Pe) 
of contained events of the and e-type. Most uncertainties cancel in this ratio and, 
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in fact, the different calculations agree to within better than 5% on the value of this 
quantity. One naively expects i? = 2 at low energies, as would result from the meson 
decay chain presented above. Detailed evaluations of R lead to a value larger than 2 
only by a few percent, for Ei, ^ 1.5 GeV. When this ratio was measured some years 
ago by the Kamiokande collaboration (Hirata et al 1988), a significant discrepancy 
with the expected value was found. This discrepancy was later confirmed by the other 
Cerenkov detector 1MB (Casper et al 1991, Becker-Szendy et al 1992b), and recently 
also by the SOUDAN II tracking calorimeter (Kafka 1993). No anomaly was found 
with the other detectors of this type, Frejus (Berger et al 1989) and NUSEX (Aglietta 
et al 1989), that are, however, statistically less significant than the Cerenkov detectors. 
These experiments report the ratio of ratios 

i?We)-|^ , (8.1) 

where Robs is the observed ratio of /i to e-type events and Rmc the Monte Carlo 
results, that take into account also experimental cuts and detector details. The values 
for this quantity are given in the table 5. 

Table 5. Ratio of observed and Monte Carlo ratios of //-type to e-type contained 
neutrino events. 



Experiment R(jt/e) = Robs/ Rmc reference 



Kamiokande 


Omtom ± 0.05^ 


(Hirata et al 1992) 




O.Q7to!^o^ ± 0.07^ 


(Fukuda ct al 1994) 


1MB 


0.54 ±0.05 ±0.12 


(Beckcr-Szcndy 1992b) 


Nusex 


o.99i°:i 


(Aglietta et al 1989) 


Frejus 


1.06 ± 0.18 ±0.15^= 
0.87 ± 0.16 ±0.08^^ 


(Berger et al 1989) 


Soudan II 


0.69 ±0.19 ±0.09 


(Kafka 1993) 



* Sub-GeV; ^ Multi-GeV Vertex contained; Fully contained. 



In figure 8 we show the momentum distribution of e-type (8a) and //-type (86) 
events in the Kamiokande data (Beier et al 1992, Hirata et al 1992), compared to a 
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Monte Carlo using the neutrino fluxes from Lee and Koh (1990). Similar results were 
also obtained by 1MB (Becker-Szendy et al 1992)). 

The deficit of //-type neutrinos (or excess of e-type, or both combined), found in 
some of the detectors, constitutes the so called atmospheric neutrino problem. If taken 
at face value, it would suggest neutrino oscillations of the type ~ Ut-, i'^ — Ve or 
v^—Vg. These would require large mixing angles (sin^2^ ^ 0.4) and, due to the typical 
V energies (GeV) and large baselines ( ^ 12000 km) they would require Am^ ^ 10~^ 
eV^. 

Another feature of the data is that there is no evidence of zenith angle dependence 
on i?(/i/e) at low energies, E^, < GeV (Hirata et al 1992), suggesting that oscillations 
should have already settled at baselines ~ 30 km for these energies. However, a 
significant angular dependence has been observed in the recent analysis of multi-GeV 
events by Kamiokande (Fukuda et al 1994) which is consistent with the sub-GeV data 
in the scenario of neutrino oscillations, because of the increase of the oscillation length 
with energy. The values of Am^ and sin^26' consistent with the Kamiokande results 
are shown in figure 9 for the — channel (9a) and i^^ — i/j- channel {9b). Also shown 
are the bounds from reactors and accelerators and from the negative Frejus result, that 
for the — Ve option exclude the parameter space allowed at present by Kamiokande. 
An explanation in terms of — oscillations would be in conflict with constraints 
from primordial nucleosynthesis on the number of effective neutrino species at T ~ 
MeV mentioned in section 5, N^, ^ 3.3 (Walker 1993), since for the required values 
of Am^ and sin^ 20 the sterile species z/g would be brought into equilibrium by the 
mixing itself (Dolgov 1981, Barbieri and Dolgov 1991, Kainulainen 1990, Enqvist et 
al 1992, Shi et al 1993) and hence lead to Ni, ~ 4. 

Since neutrino oscillations depend strongly on the neutrino energy, another way 
to test the oscillation solution to the atmospheric neutrino problem is to compare the 
upward going muons stopping inside the detector (typical energies ~ 10 GeV) 
with those going through the detector [Ej^ ~ 100 GeV). This method is independent 
of the overall normalization of the neutrino fluxes, but depends on its spectral shape. 
The analysis of the 1MB data (Becker-Szendy et al 1992a, Frati et al 1993) shows no 
evidence of an anomaly, what excludes the region with Am^ ~ 10""^ — 10~^ eV^ (since 
d ~ 10^ km) and large mixing angles shown in figure 96. 

Also the IMB-1 collaboration (Bionta et al 1988) used the non observation of an 
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anomaly between the upgoing and downgoing (E^, ~ GeV), to constrain — Vj. 
oscillations with Am^ around 10~^ eV^ (IMB-1 curve in figure 96), for which only 
the from the lower hemisphere would have oscillated sizably (and not those from 
above). 

Finally, there have been attempts to directly compare the upward-going muons 
with the theoretical expectations to further constrain the mixing parameters (these 
fluxes are consistent with no oscillations (Becker- Szendy et al 1992a, Mori et al 1991, 
Boliev et al 1991)). The uncertainties in the absolute neutrino fluxes and on the quark 
structure functions entering the inelastic cross section for muon production, make 
however those analysis less reliable (Frati et al 1993). In conclusion, the atmospheric 
neutrino anomaly could be explained in terms of oscillations between i/^ and Vr, with 
Am^ ~ 10~^ and sin^ 2d ~ 0.5. 

For the future, the results of the Soudan II experiment will significantly improve 
since their preliminary data (based on the first 1.5 kt y) still have large statistical 
uncertainties. This is important because none of the previous tracking calorimeters 
(Frejus and Nusex) observed an anomaly in the atmospheric neutrino data. Also 
relevant in this respect is the planned exposure of a Cerenkov detector to e and fx 
beams at KEK, to check for possible systematic effects due to the /i/e identification 
capability of these devices, that was never tested before. In addition, the new 
detectors Superkamiokande, SNO, Dumand, Amanda and NESTOR, that are under 
construction, will provide much more information on this issue, in particular from the 
detailed study of the zenith angle dependence of the v fluxes (Doncheski, Halzen and 
Stclzcr 1992) as well as the other methods mentioned before. As wc discussed in the 
section 6, also the proposed long-baseline oscillation experiments can help to reliably 
settle the issue of whether the atmospheric neutrinos oscillate or not. 



9. Neutrinos from Supernovae and Other Stars 

Neutrinos play an important role in stellar evolution. Neutrinos produced in 
the interior of stars can just stream out, except at the large densities of collapsing 
supernovae, where they are trapped for some time. While in ordinary main-sequence 
stars, such as our sun, the photon luminosity is larger than the neutrino luminosity, 
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for stars whose central temperature is higher than 10 keV (O and Si burning stars 
and beyond) neutrino emission is the main form of energy foss. Thus, stars can test 
any non-standard neutrino property that increases the neutrino energy foss rate from 
stars or the neutrino energy transfer inside stars, thus changing their standard course 
of evolution (for a review see Raffelt 1990a). For example, an increase of energy loss 
accelerates the evolution of stars, shortening their lifetimes. 

Nowhere the effect of neutrino cooling is more dramatic than in a supernova 
explosion. The observation of the neutrino burst from the supernova SN1987A 
provided a confirmation of the theoretical understanding of the stages of the 
gravitational collapse of a type II (core collapse) supernova (see e.g. Burrows 1990 
or Hillebrandt and Hoflich 1989). Most of the energy liberated in the collapse of the 
iron core of the progenitor star into a neutron star, i.e. ~ 3 x 10^^ ergs, is emitted in 
neutrinos in about 10 sec. While the collapse is well understood, the actual explosion 
which ejects the outer layers of the original star is not. Because only ~ 1% of the total 
energy liberated in the collapse goes into the explosion, relatively very small eff'ects 
have to be well understood to account for it theoretically. Therefore, supernovae 
explosions are very difficult to obtain in simulations. This difficulty is sometimes 
called the "supernova problem" . Neutrino properties that could transfer some of the 
neutrino energy to the outer layers of the star have been repeatedly advocated as a 
solution to this problem. 

A surprising wealth of information on neutrinos have been obtained from SN1987A, 
considering that only 19 neutrino events were observed (by two large water Cerenkov 
detectors, 1MB, Bionta et al 1987 and Bratton et al 1988, and Kamiokande II, Hirata 
et al 1987 and 1988, through the reaction u^+p ^ n + e"*"). Let us mention some of 
the bounds obtained (for a review see e.g. Raffelt 1990a and 1990b). 

The simple arrival of (electron anti-) neutrinos from the SN1987A, provides a 
bound on any effect that could have removed them from the burst before arriving 
to earth, such as decay, scattering, or interaction with magnetic fields. The lower 
bound obtained on their lifetime is Ti^^ > 6 x 10^ sec (m^^/eV). Due to the possibility 
of a difference between current and mass eigenstate neutrinos this bound does not 
exclude that part of the electron neutrinos emitted in the sun could decay before 
arriving to the earth (Prieman, Haber and Preese 1988), explaining in this way the 
solar neutrino problem. However decaying neutrinos do not provide a satisfactory fit 
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to the solar neutrino observations (Acker and Pakvasa 1994). The neutrinos arriving 
from SN1987A where not removed by scattering, for example, with a background 
of relic Majorons, thus the neutrino majoron coupling (see section 3) must be 
{9ee)u^u^j < 10~^ (Kolb and Turner 1987). The arrival of neutrinos provides also 
a bound on their charge. If neutrinos would have a non-zero electric charge Qiy^, their 
trajectories would have been bent in the galactic magnetic field. The observation of 
the neutrino burst gives a limit Q,^^ ^ 3 x 10~^^e (Barbiellini and Cocconi 1987) (while 
the best limit from the neutrality of atoms is Qi^^ ^ 10~^^e, Bauman et al 1989). 

An upper bound on rrii,^ is due to the (non) dispersion of the neutrino emission 
time. The neutrino signal was spread over ~ 10 seconds and the emission time could 
not have been much longer. Massive neutrinos from the star arrived to earth with a 
delay time (time of arrival minus time of emission) At = d/v ~ d[l + {l/2){m^/ Ejy)'^] 
where d ~ 50 kpc= 1.6 x 10^ lyr is the distance travelled and E^, is the neutrino 
energy. If all neutrinos (with equal mass m^) had been emitted at the same time, the 
more energetic ones should have arrived earlier. This is not what happened (energies 
and times were not so correlated). Thus, by increasing the assumed mass rrii,, the 
inferred dispersion of the emission times increases and for some mass value becomes 
unacceptably large. This argument yields the upper bound rrij^^ < 23 eV (at 95% C.L. 
ignoring systematic uncertainties, Loredo and Lamb 1989). This bound (as those of 
section 4.1) is independent of the Dirac or Majorana nature of the neutrinos. There 
arc no bounds on the other flavour neutrino masses because the bulk of the events 
observed were z^g, due to detection cross-section arguments. 

The number and energy of the detected events correspond well with what is 
expected, leaving not much room for new sources of energy loss. Also a new mechanism 
for energy loss would have speed up the cooling, thus shortening the pulse with respect 
to what was seen. For example Majoron models (singlet or mostly singlet Majorons, 
the only viable ones, see section 3) are constrained by these argument, that reject a 
region in rxij^-V space, for V , the scale of spontaneous symmetry breaking, in the range 
few GeV ^ 1 TeV (Choi and Santamaria 1990). If neutrinos are Dirac particles, 
their mostly sterile right-handed helicity components can be produced due to their 
m/E admixture in the left-handed chirality component (the interactive component). 
Subsequently they escape, because their cross sections are too small to trap them, 
increasing the cooling rate of the neutron star. A bound < 0(10 keV) has been 
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obtained (RafFelt and Seckel 1988, Burrows, Gandhi and Turner 1992, Mayle et al 1993 
and references therein) from preventing a shortening of the duration of the observed 
coohng neutrino pulse. This bound apphes only to neutrinos lighter than 1 MeV, that 
can be copiously produced, because the production of heavier neutrinos within the star 
is approximately suppressed by a Boltzman factor (Sigl and Turner 1994). The upper 
bound of 0(10 keV) is hard to improve due to the difficulty in evaluating the effect of 
the very fast multiple scattering of nucleons inside the supernova on the emission of 
the right-handed neutrinos (Raffelt and Seckel 1991 and 1992). This difficulty affects 
most bounds (also on other particles) based on energy loss in SN1987A. Sterile right 
handed neutrinos could also have been produced due to new interactions with an 
effective coupling constant Gr (the equivalent of the usual Fermi constant Gp)- The 
above cooling arguments yield Gr ^ IQ-^Gf (Raffelt and Seckel 1988, 1991, 1992). 

If neutrinos are Dirac particles and they arc unstable, the that escaped from 
the inner, hotter, regions of the supernova, could decay into energetic [E ~ 100 
MeV) left-handed neutrinos that would have been seen as extra events if 1 keV 
^ mix ^ 300 keV and 10~^sec(m/l keV) ^ r ^ 5 x 10^sec(m/l keV), thus these ranges 
are forbidden for Dirac neutrinos (Dodelson, Frieman and Turner 1992). Following a 
similar argument, one finds a more dubious bound that would apply to both Dirac 
and Majorana neutrinos in the same mass range, emitted from the "neutrino sphere" , 
the boundary of the region where neutrinos are temporarily trapped (so these are 
interacting neutrinos), if electron neutrinos are produced in the decay. Then, lifetimes 
3 X 10^ sec (keV/m) ^ r ^ 2 x lO^'^ sec (keV/m) are rejected (Gelmini, Nussinov 
and Pcccci 1992, Mohapatra and Nussinov 1992, Soares and Wolfenstein 1989). This 
bound is less meaningful because neutrinos are emitted from the neutrino sphere with 
energies of 0(10 MeV), consequently the energy of the decay products is smaller than 
before, thus reducing the expected number of events and their energies, what makes 
these events more difficult to distinguish from background events. These bounds on 
lifetimes apply to invisible decay modes, in which the decay products do not include 
photons or e+e~, i.e. v — > Su' or — > u'cj), where is a Goldstone boson, a Majoron in 
most models (see section 3). Much more restrictive bounds apply to radiative decays, 
as shown below. 

Turning now to the electromagnetic properties of neutrinos, a magnetic or electric 
dipole moment fii, would also allow for helicity ffips into otherwise inert vr of Dirac 
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neutrinos in the electromagnetic scattering of the trapped vl with charged particles. 
Thus, the same arguments on energy loss from SN1987A impose, for Dirac neutrinos 
only, < 10~^^iJ,B, with the uncertainty mentioned above (Barbieri and Mohapatra 
1988, Lattimer and Cooperstein 1988, Goldman et al 1988) for neutrinos that can 
be copiously produced in the supernova, namely m ^ 1 MeV. Here Hb = e/2me 
is the Bohr magneton. Other bounds, valid for nii, ^ 10 keV, arise from energy 
loss due to the plasmon decay 7 pffj in red-giants and in white dwarfs. These 
are respectively jij, < Z x 10"^^/^^ (Raffelt 1990c, Raffelt and Weiss 1992) and 
111, < 1 X 10~^Vb (Blinnikov 1988, Raffelt 1990a). These two bounds apply to both 
Dirac (direct or diagonal and transition) and Majorana (transition) neutrino dipole 
moments ('diagonal' and 'transition' refer to flavour space, see section 2). Note that 
these bounds would exclude the spin-flip solution to the solar neutrino problem (see 
section 7). 

For comparison, let us mention here the cosmological bound based on the 
effective number of relativistic neutrinos during nucleosynthesis, Ni, < 3.4 (Olive 
et al 1990, see section 5). It requires /ii, ^ 5 x 10~^^ fiB for Dirac neutrinos 
(Morgan 1981, Fukugita and Yazaki 1987), because otherwise the mostly inert right- 
handed helicity components could be brought in equilibrium in the early universe 
due to the electromagnetic interactions of left-handed neutrinos. On the other 
hand, the laboratory bounds for magnetic moments (Review of Particle Properties 
1994 and references therein) are n^,^ < 1.08 x 10~^/Ub, h^,^ < 7.4 x 10~^^fXB and 
< 5.4 X 10~^hb- This latter bound on ii^,^ applies only to a direct moment 
(Coopcr-Sarkar et al 1992). Using the same data ( Cooper- Sarkar et al 1985), 
Babu, Gould and Rothstein (1994) found a laboratory bound for transition moments 
A*z^T-tran < 1-1 X 10~^ (McV/mj^^ )//b 5 that is more restrictive for an MeV tau neutrino. 
In the simplest extension of the SM, where only uji and Dirac masses m,y are added 
(section 3.2), the predicted neutrino magnetic moment is much smaller than any of 
these bounds, //j^ ~ 3 x 10~^^ {mi,/eV)iJ,B (Lee and Schrock 1977, Fujikawa and Shrock 
1980). 

The non-observation (by the Solar Maximum Mission Satellite, SMM) of a 7- 
ray burst in association with the SN1987A neutrino burst, forbids radiative neutrino 
decays u — > ^u' with lifetimes Trad > 2 x lO"*^^ sec {m,^/eV) for m,^ < 20 eV. For 
larger masses, the photons produced would be spread out in time and the upper bound 
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weakens to Trad ^3x10 ^ sec for 20 eV ^ m„ ^ 100 eV and Trad > 0.8 x 10^^{eV/m„) 
for 100 eV ^ m^, ^ 1 MeV (Kolb and Turner 1989, Raffelt 1990a and 1990b 
and references therein). The mass-hfetime region excluded by the SMM bounds is 
shown in the hatched area of figure 2. These bounds on Trad translate into very 
restrictive bounds on transition moments, for all neutrino types and flavours (through 
the relation Trad = 87r//ij^^m^^, valid for negligible flnal neutrino mass). These bounds 
are |U^tran < 1 X 10~^(eV/mj^)^|UB for ^ 20 eV and //i,tran < 5xl0~^^(l eV/mj^)iiB 
for 20 eV ^ <, 1 MeV (RaflPelt 1990a). 

The mass vs. lifetime bounds for visible decay modes of neutrinos heavier than 
1 MeV (only possible for v^.) have been recently reanalysed (Sigl and Turner 1994, 
Babu, Gould and Rothstein 1994). These neutrinos have another visible decay mode 
besides v^- v'^, namely v^- — > v'e^e~ . Even if this is the dominant decay, the SMM 
bound can be applied, because photons will be produced with branching ratio 10~^ 
through bremsstrahlung (Sigl and Turner 1994) or by other processes (Mohapatra, 
Nussinov and Zhang 1994). The SMM bound only applies if neutrinos decay outside 
the supernova, i.e. for r > 100 sec. For shorter lifetimes, the decay products would be 
trapped inside and bounds result from their effect on the supernova energetics. The 
combinations of both bounds exclude a massive v^- of lifetime 10~^ sec ^ t ^ 10^ 
sec decaying into visible modes (Sigl and Turner 1994). The excluded region for 
rrij^ > 1 MeV taken from Sigl and Turner (1994) is included in the hatched area of 
figure 2. Laboratory bounds (Cooper-Sarkar et al 1985, Babu, Gould and Rothstein 
1994) exclude these modes for even shorter lifetimes, r ^ 0.1(m;y/MeV) sec *. So, 
if the dominant decay modes are visible, all lifetimes shorter than ~ 10^ sec are 
forbidden by the combination of SN1987 bounds and laboratory bounds. If we add the 
nucleosynthesis bounds (see section 5, Dodelson, Gyuk and Turner 1994, Kawasaki et 
al 1994) that forbid Ur heavier than ~ 0.3 MeV if the lifetime is shorter than ~ 100 sec, 
the combination of the three bounds reject a Vr heavier than 0.3 MeV that decays 
dominantly into visible modes for any lifetime. 

The MSW effect (see section 7) could also happen within a supernova. The 
oscillation of or 1/^ into i^e outside the neutrino sphere could have important 
consequences. The reason is that the emitted and Vr are more energetic than 

* Very unlikely short lifetimes t ^ 2 x 10~^^(m^/MeV) would be allowed only for the mode Vt — > t's7) 
where i/s is a sterile neutrino. 
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the emitted v^- This is because and only interact through neutral currents with 
the surrounding matter. This means, they have a larger mean free path than (that 
also interact through charged currents) and they consequently emerge from deeper 
and, thus, hotter layers of the supernova core. Thus, if a large fraction of the Vr (or 
z/^) are transformed into energetic (that, as we just said, have larger interaction 
rates with the surrounding matter), the increased energy given to outer layers of the 
star may help explode them (Fuller et al 1992). However, the same flavour conversion 
may preclude nucleosynthesis through r-processes in supernovae. Heavy elements are 
believed to be synthesized through rapid neutron capture processes (r-processes) in 
the neutron-rich outer layers of an exploding supernova. The more energetic electron- 
neutrinos, resulting from a MSW conversion of outgoing zv^ or v^- before reaching the 
region where r-processes should occur, would reduce the amount of neutrons in the 
material of the crucial layers (through i/g + n — p + e~). The environment would then 
become proton-rich and this would stop r-processes entirely. Because there is no other 
known production site for the elements that should be produced though r-processes in 
supernovae, this effect must be avoided. This excludes the region of the Am^ — sin^ 20, 
for z/e — and z/e — i't oscillations, corresponding to approximately Am^ ^ 4 eV^ for 
sin^ 29 ^ 10~^ (Qian et al 1993). This is an important region, because it corresponds 
to a Vr or with a mass suited to be hot dark matter (between a few eV and a few 
tens of eV, see section 5), if v^. is much lighter. 

The 1MB and Kamioka cerenkov detectors have shown the feasibility of observing 
the neutrinos from supernova explosions, and as we just discussed their few events 
allowed to obtain significant information on neutrino properties. In the future, with 
the construction of Superkamiokande and also especially of the SNO detector sensitive 
to neutral currents, the detection of a supernova within our galaxy (there are ~ 2 
galactic supernova explosions per century!) may allow to constrain the masses of 
z/^ and v^- if they are larger than ~ 25 eV (Minakata and Nunokawa 1990, Seckel, 
Steigman and Walker 1991, Krauss et al 1993, Burrows, Klein and Gandhi 1992), and 
maybe even down to 15 eV with new type of detectors (Cline et al 1994). It may 
also become possible to detect the background of neutrinos from all past supernovae, 
that should be the dominant neutrino flux at energies just beyond those of the solar 
neutrinos (20-30 MeV). 

Finally, let us mention the under-water (ice) neutrino telescopes under construction 
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at present with energy thresholds between a few and 10 GeV. They will observe the 
high-energy neutrinos that should be emitted in comparable numbers with 7 rays in 
the most energetic galactic and extragalactic processes, such as the shock acceleration 
in Active Galactic Nuclei. These experiments, DUMAND, AMANDA, Baikal and 
NESTOR, with an area of approximately 0.02 km^, can be considered the prototypes 
for the final desirable full size experiments of 1 km^ (for a review see Halzen 1993). 



10. Concluding Remarks 

The hints for non-zero neutrino masses in solar and atmospheric neutrinos, and 
the possibility of confirming or rejecting them in the near future, make the subject of 
neutrino masses a particularly exciting field of research at present. So much so that 
one could indulge in the actually premature exploration of the consequences of the 
confirmation of both. One could even wonder about their compatability with other 
possible desirable mass values, such as masses of order 1-10 eV for neutrinos to account 
for part of the Dark Matter in the Universe, or Majorana masses within the reach of 
present neutrinoless double beta decay experiments, i.e. < m^^ >~ 1 eV. 

The MSW solution to the solar neutrino problem would require ml. — m^^ ~ 
O(10~^) eV^, while the solution to the atmospheric neutrino deficit requires — 
mf,^ ~ O(10~^) eV^. An obvious simultaneous solution for both requires mj^^ ~ 
10~^ eV, mjy^ ~ 10~^ eV and rrii^^ -C rn^^, with Ui = and Vj = Ut. These mass 
values can be obtained in a see-saw model with a large right-neutrino Majorana mass 
M ~ 10^^ GeV (see section 3). 

Incorporating either the solar neutrino solution or the atmospheric neutrino 
solution with either DM neutrinos or {m,^J ~ 1 eV is easy in see-saw models with 
smaller values of M (see e.g. Peltoniemi, Tommasini and Valle 1993). A solution 
incorporating both the solar and the atmospheric neutrino solutions and either DM 
neutrinos or {rriv^) — 1 eV requires quasi-degenerate neutrinos, where the mass 
differences necessary to account for the oscillations are much smaller than the masses 
themselves (Caldwell and Mohapatra 1993, Peltoniemi and Valle 1993). Models of 
this type, possibly with extra inert neutrinos, can even incorporate all four mass 
requirements mentioned. In a possible solution the three active neutrinos are almost 
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degenerate. This can be obtained if the see-saw mechanism determines the mass 
differences but the left-neutrino Majorana masses give the masses themselves (see 
section 3). 

The confirmation of a see-saw mechanism is by no means the only possible outcome 
of neutrino mass experiments in the near future. We may uncover the Majorana nature 
of neutrinos in future neutrinoless double beta decay experiments if the effective 
mass is not much smaller than 1 eV (section 4.2). Still, we may find neutrino masses 
in the range forbidden by the present energy density of the universe (in direct mass 
searches on with a galactic supernova explosion), what would be a strong indication 
of the existence of Goldstone bosons (typically singlet or mostly singlet Majorons, see 
section 3), necessary for the neutrinos to decay or annihilate fast in the early universe 
(section 5). Or even neutrinoless double beta decay with the emission of a boson, that 
would be an even more direct indication of Goldstone bosons. Even if these results 
are possible, they are unlikely. We do expect instead a resolution of the solar and 
atmospheric neutrino problems within the next decade. These are certainly exciting 
times for neutrino physics. 
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Figure captions 



Figure 1. Bounds on unstable neutrinos that apply on any decay mode. The gray 
area is rejected. The continuous contour shows the bound flop < ^o- The dashed 
contour shows the structure formation bound only for experimentally allowed masses. 
The dot-dashed contour shows the nucleosynthesis bound on the allowed number of 
effective extra neutrino species 5Ni, (see section 5). 

Figure 2. Bounds on unstable neutrinos whose main decay products include photons 
or e"'"e~ pairs. The gray area, excluded by cosmological limits (distortions of the CBR 
and other photon backgrounds), is reproduced from the figure 5.6 of Kolb and Turner 
1990. The hatched area is excluded by the non observation of photons in coincidence 
with the observed neutrino flux from the supernova SN1987A (see section 9). 

Figure 3. Present constraints on Am^ — sin^ 26 from reactor and accelerator searches 
of — Un oscillations (a) and for j/^ — u^- (b). For references see the text. 

Figure 4. Sensitivity of future oscillation searches together with the present bounds 
and the regions suggested by atmospheric neutrinos and to account for the dark matter. 
Figure 4a is for — oscillations while 46 is for v^^ — v^-. 

Figure 5. Standard solar neutrino spectrum (Bahcall and Ulrich 1988) and 
experimental thresholds for the existing detectors. 

Figure 6. Survival neutrino probabilities vs. E/l^m? for sin 2^ = 0.1. 

Figure 7. MSW contours for the different experiments and preferred values of 
Am^ — sin^ 29 (shadowed). Also shown is the region excluded by the non-observation 
of a day-night asymmetry by Kamiokande. 

Figure 8. Spectrum of e-type and /U-type atmospheric neutrinos measured by 
Kamiokande. The histogram is the MC prediction using the fluxes of Lee and Koh. 

Figure 9. Values of Am^ — sin^ 29 suggested by Kamiokande measurement of R{jj,/e) 
atmospheric ratio and its zenith angle dependence, together with the bounds from 
accelerators and reactors for i/^ — Ue (a) and — Vr (b) oscillations. Also shown are 
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the regions excluded by other searches with atmospheric neutrinos: Frejus value of R, 
IMB-1 up/down i/^, and 1MB stopping/through-going muons. 
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